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ABSTRACT 



An experimental investigation was conducted to determine 
the dynamic responses of stiffened flat plates to impulsive 
loading by underwater shock waves . Air-backed flat plates 
with machined external stiffeners and clamped boundary 
conditions were subjected to shock loadings from TNT charges 
detonated underwater. The plates were instrumented to 
measure transient strains and free field pressures. Test 
results were compared to preshot calculations done using the 
USA-STAGS code. Particular emphasis was placed upon the 
code 's ability to predict stiffener tripping and shear at 
plate boundaries. 
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INTRODUCTION 



I . 



A. BACKGROUND AND HISTORY 

Sergei Gorshkov, Admiral of the Soviet Fleet, has 
observed that the end of World War Two marked, "... 
the start of a military-technical revolution which in scope 
and depth transcended all the reforms and transmissions 
which had previously occurred in the armies and fleets of 
the world." [Ref. 1] 

Much has been made of the tactical and geo-po 1 it ical 
implications of that revolution. Less popularly publicized, 
but equally radical are the changes that the developments in 
technology have spawned in the processes of ship design and 
const ruct ion . 

Throughout most of history , war at sea has been a 
process of trying to set one's opponent afire, board him, or 
hurl chunks of stone or metal at him across relatively short 
distances. All of these methods had the distinguishing 
characteristic of being essentially disabling tactics, 
unless a ship was burned or perforated to the point of 
sinking (and they often were not), damage was generally 
largely superficial and the ship could be refitted to see 
service again. Witness the many vessels that have been 
fought on both sides of a conflict and the tradition of 
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expecting naval personnel to supplement their pay with prize 
money . 

Two events occurred at or around the beginning of the 
twentieth century to change all of that: the invention of 

high explosives and the advent of the submarine. 

Initially, the submarine was primarily a commerce 
raider. It was too small to carry a prize crew and too 
lightly armed to effectively combat a conventional warship. 
During World War One, most battles between submarines and 
conventional ships were fought on the surface. Submarines 
would use their stealth to sneak-up upon unsuspecting mer- 
chantmen, surface and then attempt to sink them with their 
deck cannons. Torpedoes were largely unreliable and the 
submarines themselves were at risk if fast armed vessels 
were present. Convoying the merchants with small, agile 
men-of-war was, briefly, an effective means of countering 
the submarine threat. However, as torpedoes became more 
effective so, too, did the submarine. For the convoy escort 
vessels it was no longer sufficient to have to wait for a 
submarine to surface in order to attack it. An effective 
standoff weapon was needed. 

That standoff weapon was the depth charge. First intro- 
duced in World War One, the depth charge came of age in the 
Second World War; and with it a new era in ship design was 
ushered into existence. For the principle of this new 
weapon was not one of punching a hole in a ship's hull and 
then exploding internally; rather, high explosives permitted 
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generation of intense pressure waves in the sea that could 
break a vessel's back or rupture hull fittings whether the 
vessel's skin was punctured or not. Simply increasing the 
thickness of plating in discrete areas of vulnerability was 
no longer sufficient. Design for strength and resiliency 
under unequal pressure loadings was required. Moreover,' 
this requirement applied to attacker as well as attacked. 
Although the attacking surface vessel had the advantage of 
being able to run from the immediate scene of the explosion 
and the effects were not compounded by as great a pre- 
existing static head of water, these weapons were not direc- 
tional and could be expected to damage all within their 
sphere . 

The issue of a vessel's ability to withstand an under- 
water explosion, be it submarine or surface ship, became 
more pressing as the explosive devices became increasingly 
powerful and sophisticated. Not only did depth charges 
become lethal over a wider area, but the advancing tech- 
nology permitted mines, too, to grow from essentially 
contact weapons to devices that could endanger all types of 
shipping from considerable depths. The issue has become 
particularly crucial with the development of nuclear weapons 
which are so powerful that no device need contact a ship or 
group of ships to cause rampant destruction. 

The response of the U.S. Navy was to require that ships 
and their critical components be designed and built to 
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withstand specified levels of shock loading. These re- 
quirements included (and continue to include) actual shock 
testing of the first ship of a class and many crucial systems 
on the component level. 

During the immediate post-World War Two period, many of 
these requirements were largely tentative and experimental. 

By the early 1960's, however, sufficient experience had been 
gained to allow design requirements to be set forth 
[Ref. 2]. These specifications established specific guide- 
lines for contractors to meet based upon type of equipment, 
type of vessel in which the equipment was to be used, and 
location and method of mounting within the vessel. As with 
the ships themselves, actual shock testing of the item was 
the primary means of establishing compliance; as it remains 
today . 

Performing a shock test upon production end products-, 
particularly when' those items are as large as an aircraft 
carrier can have some obvious drawbacks: it is costly, it 

is dangerous, and the results can be embarrassing. 

Moreover, there can be longterm structural damage which 
is neither immediately apparent nor correctable, but which 
does not by itself constitute non-compliance of contract 
terms. Moreover, some items simply do not lend themselves 
to shock testing. 
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What has been needed is a method of analyzing the design 
itself in order to achieve a numerical approximation of the 
end product's response to shock. 

This need has coalesced into the development of two 
computer codes; EPSA (Elasto Plastic Shell Analysis), 
produced by Weidlinger Associates, and USA-STAGS (Underwater 
Shock Analysis--Structural Analysis of General Shells) 
produced by Lockheed Missiles and Spacecraft, Inc. 

B. OBJECTIVES 

This study is the second in an ongoing series sponsored 
by the Defense Nuclear Agency (DNA) into the effects of 
underwater explosions upon ship hulls. The intent of these 
studies has been to test the applicability of the two previ- 
ously mentioned structural response, underwater shock anal- 
ysis computer codes to naval design and analysis requirements. 
Under the umbrella of general applicability, focus has 
been placed upon two defined goals: 

1. Use of the codes to gain insights into large deflec- 
tion elastic-plastic responses of a submerged struc- 
ture subjected to transient acoustic shock loading; 
with especial emphasis upon stiffener tripping 
phenomena . 

2. Performance of underwater explosion testing of the 
structure to validate the results predicted by the 
code . 

To achieve these goals, flat plate models geometrically 
similar to a ship's stiffened hull structure have been 
studied. These plates were air backed to simulate general 
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conditions that could be expected in a ship's plating; and 
were subjected to shock loadings that produced deformations 
well into the plastic regime in order to ensure stiffener 
tripping. Existing information concerning underwater explo 
sion phenomena was employed in an attempt to reduce experi- 
mental uncertainties. 

The first study in this series was that conducted by LT 
Thomas Rentz, USN in fulfillment of thesis requirements at 
the Naval Postgraduate School [Ref. 3] and [Ref. 4]. In 
that study he developed the following experimental strategy 
to study the EPSA code: 

An attempt to validate the code was conducted in two 
phases: Using pressure approximations based upon empiri- 

cally determined formulae, the code was first used as a 
pre-test predictor of test results. Once the underwater 
explosion testing had been conducted, pressure histories 
derived from that testing were input into the code. The 
plate strains predicted by the code in both its preshot and 
post-shot capacities were then compared to actual strain 
histories from testing. Plate deformations and boundary 
stresses predicted by the code were compared to observed 
physical responses. 

In this study, the code to be tested is USA-STAGS. In 
an effort to maintain continuity with the earlier study in 
the series, the same experimental strategy will be followed 
The test geometry and flat plate model will also be 
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retained. Conclusions regarding the ability of USA-STAGS to 
predict the responses of a simple model to underwater explo- 
sions will be drawn; and, perhaps of even greater impor- 
tance, guidelines will be established for using the 
USA-STAGS code and conducting underwater explosion testing 
for future studies. 
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II. UNDERWATER EXPLOSION THEORY 



A. THE PRESSURE PULSE 

Any conversion of matter, whether chemical or nuclear, 
which results in the very rapid production of large quanti- 
ties of gas at very high temperature can be categorized 
within the broad species "explosion." Typically, for mili- 
tary high explosives such as TNT the pressure and tempera- 
ture within the gases produced by an explosion are on the 
order of 50,000 atmospheres and 3000 degrees Celsius 
[Ref. 5: p. 3]. For nuclear explosions the initial tempera- 

ture is on the order of a million degrees Celsius and the 
pressure, which is dependent upon the yield is similarly 
orders of magnitude higher than that for a conventional 
explosive [Ref. 6]. 

Clearly, temperatures and pressures of the levels 
described must dissipate high levels of energy through any 
surrounding medium. In the case of water, the initial mani- 
festation of this disturbance is an intense compression 
created pressure wave propagating radially outward from the 
charge; followed nearly instantaneously by a displacement of 
the fluid. In the immediate vicinity of the charge, the 
velocity of propagation of the pressure wave is several 
times the speed of sound (about 5000 feet per second in 
water) but approaches the value almost immediately and 
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therefore may be apprehended as an essentially acoustic 
phenomenon. As the water begins to be displaced, the gas 
bubble created by the explosion expands; thereby decreasing 
its internal pressure and correspondingly the pressure on 
the constraining fluid. Viewed on a pressure-time curve the 
initial rise in pressure is so nearly vertical as to appear 



discontinuous, while the pressure relief is characterized by 
an approximately exponential decay (Fig. 2.1) [Ref. 5: pp . 

4-7] . 




Figure 2.1 Simple Exponential Approximation of Incident 
Shock Wave 



The actual magnitudes of the pressure pulses and the rate 
of exponential decay have been found through experiment to be 
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functions of the type of explosive, the charge weight, and 



decay constant were originally determined by Robert H. 

Cole and were presented in his work, Underwater Explosions 
[Ref. 5]. These have since been simplified to use a series 
of empirically derived coefficients developed by R. S. 

Price [ Ref . 7 ] . 

The pressure profile can be expressed by the equation 



where P q is the initial (greatest) pressure of the shock 
wave, 9 is the decay constant describing the exponential 
decay, and t - t^ is the time elapsed from the arrival of 
the shock front . 

The initial pressure, P , and the decay constant, 

9, are dependent largely upon the type of explosive and 
the weight of the charge. These may be expressed 



the standoff distance.' 1 ' The equations for the pressure and 



P ( t ) = P exp 
o ^ 



( eqn 2.1) 




( eqn 2.2) 



and , 



9 




( eqn 2.3) 



''"The standoff distance is defined as the distance from 
the charge to the nearest point on the target. 
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where P q is in pounds per square inch and 3 in milli- 
seconds. K 1 , K- 2 , , A 9 are the explosive dependent values 

(Table I). W and R are the charge weight (in pounds) and 
standoff distance (in feet), respectively. 

Table I 

Explosive Constants I 



Explosive 


TNT 


HBX-1 


PENTOLITE 


K x 


22505 


22347.6 


24589 


N 


1 . 18 


1 . 144 


1 .194 


K 2 


.058 


.056 


.052 


A 2 


- .185 


-.247 


- .257 



B. SURFACE CUTOFF AND BOTTOM REFLECTION 

Shock waves propagate radially from the charge. Thus, 
for a point source in an infinite medium, the wave could be 
expected to travel as an ever expanding sphere until damped 
by the fluid. For small to moderately sized charges deto- 
nated well below the surface in great depths of water, this 
is a satisfactory analogy. For explosions nearer the 
surface or in shoal water, however, there are secondary 
effects for which there must be an accounting. 

Y(hen a shock wave ’reaches the surface of the water, the 
adjoining air is of sufficiently lesser density to be inca- 
pable of supporting the shock wave. The result is that the 
effective pressure at the interface must be zero in 
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This implies that 



compliance with the laws of continuity, 
a negat ive shock wave of magnitude equal to the original 
incident wave must be propagated downward through the 
medium to satisfy equilibrium conditions. This negative 
wave will travel through the water at an angle to the 
surface that is equal to that of the incident wave; much 
as the negative reflection image from a mirror continues at 
an angle equal to the angle of incidence. 

Remembering that these waves are radiating spherically 
outward from the point of inception, it is clear that both 
the incident wave and the generated surface wave will even- 
tually reach the ship, submarine, or other target object. 

It is likewise clear that, with its longer path, the surface 
wave must reach the target at some time later than the inci- 
dent; which travels only the shortest, or standoff, distance 
between charge and target. The net effect is additive. 

When the surface wave reaches the target, the initial shock 
front will have passed and the incident wave will be at some 
point in its exponential decay. The negative surface wave 
will then "chop off" the tail of the incident wave. If the 
spatial geometry of the shot is such that the magnitude of 
the surface wave is greater than the remaining magnitude of 
the incident wave, a region of negative pressure will be 
created (Fig. 2.2). If this negative is less than 
the pressure required to keep the fluid in a liquid state, 
the fluid will "flash" into vapor. This chopping effect is 
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referred to as surface cutoff, and the negative pressure as 
bulk cavitation. 




The direct effect of bulk cavitation is, by itself, not 
particularly harmful. The structure simply experiences 
essentially no pressure while it is in the cavitated region. 
When, however, the static head of water above the region 
overcomes the cavitation forces, the region closes suddenly 
to generate a destructive reloading upon the structure. 

To further complicate matters, if the shot occurs in 
shoal waters or at depths near the bottom a second boundary 
is also present. Unlike the surface, in this case the 
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incident wave is impinging upon a medium of greater density 
than the fluid. The reflected wave will therefore be 
positive. The actual magnitude of this wave depends upon 
the physical characteristics of the bottom material. The 
effects of this bottom reflected wave are also additive 
and will generally reach the target at some time after 
the initial shock front has passed; although the actual 
time of arrival of this secondary front will once again 
be dependent upon the bottom conditions and the speed at 
which the wave travels through the bottom medium. Addi- 
tionally, in shoal waters, that same bottom reflected wave 
will reach the surface; creating its own surface cutoff. 

The net effect upon the target is thus one of an inci- 
dent shock wave; which, if the charge and target are suffi- 
ciently near the surface, will be cutoff by the negative 
surface wave. This, in turn, is followed by the bottom 
reflected wave which also experiences an exponential decay, 
and which may also be cutoff by its own generated surface 
wave ( Fig . 2.3). 

Analytically, the easiest way to model the effects of 
incident, surface, and bottom waves upon the target is as 
has been represented in Figure 2.4, with three separate 
images . 

The incident and surface images are relatively easy to 
depict. The incident image, iVL , of course represents the 
charge itself and occupies the same location in the spatial 
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H o-_ 



SURFACE 




Figure 2.4 Three Image Presentation of Shock Waves 
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the same magnitude 



geometry. The surface image, M , is of 
as the incident image but of an opposite sign. It lies on a 
line drawn vertically through the charge at a distance above 
the surface equal to the depth of the charge below. 

The surface cutoff pressure may therefore be 'computed 



p = _ K |™ll 3 V Al 

c 1 R 



where R is now the distance 
s 

the target. The time of cuto 
initial shock front i's thus t 



to cover the greater distance 



t = 
c 



R - R. 
s 1 



( eqn 2.4) 

from the surface image to 
ff after the passage of the 
he time difference for the wave 
, or 

( eqn 2.5) 



where c is the speed of sound in the fluid. The cutoff 
pressure will then decay in the same exponential manner as 
the incident pressure. 

As may be conceived from the preceding discussion of the 
bottom reflection; the location of the bottom image is not 
so straightforward . It may be viewed as initiating along 
the same vertical line as do the incident and surface 
images, and at a depth below the fluid-bottom interface 
equal to the distance of the charge above. Here, however, 
the similarities to the surface image end. The bottom image 
cannot be considered a point charge because of the tendency 
of the pressure wave to "smear" across the bottom. The 
magnitude and speed of the wave, too, are dependent upon 
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bottom characteristics. A muddy bottom, for instance, will 
tend to reflect a diffuse, low magnitude wave travelling at 
or about the speed of sound in water. A rock bottom on the 
other hand will produce sharp, high magnitude pulses which 
will tend to radiate over a distance due 'to the higher speed 
of travel of the shock wave in the bottom medium than in the 
fluid.' Due to this dependence upon bottom conditions, 
determination of bottom reflection effects do not lend them- 
selves to general predictive studies and are often disre- 
garded, although they can be significant. 

C. BUBBLE EFFECT 

The foregoing ' discussion of the development of the pres- 
sure wave and its impact upon the target assumes that a 
single, clean wavefront is produced by the explosion. In 
point of fact, such is rarely the case. When the charge 
explodes, it first creates a high intensity shock front 
which is the initial incident shock wave. It also produces 
a bubble of very hot gasses at very high pressure. 

The bubble expands until it reaches pressure equilibrium 
with the surrounding water. When it does so, a portion of 
the energy within the bubble is released as a new shock 
front . This release of energy then permits the momentum of 
the water displaced to assert itself and the bubble is 
collapsed, compressing the gases within. The displaced 
water can be conceived as attempting to regain its original 
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pre-bubble condition. As the water rushes in, it compresses 
the gases beyond the equilibrium condition and the bubble 
again expands, again releasing a portiop of its energy as a 
shock wave, albeit a smaller one than the preceding (Fig. 
2.5). This process has been known to produce as many as ten 
significant pulses [Ref. 5: pp . 8, 9], which, though weaker 

than the initial shock front, can still have a pronounced 
effect upon a structure which has already been loaded by the 




Figure 2.5 Multiple Expansions of the Hot Gas Bubble 
[Ref. 5 : p . 7 ] 



29 



The time and radius of the bubble at the first pulse 
has, again, been determined empirically, the maximum radius 
as 



A = K„, 
max 3' 



W 



1/3' 



S D + 33) 



1 / 3i 



( eqn 2.6) 



at time, 



w 1/3 _ ^ 

(D + 33 ) 5/6 J 

where A is expressed in feet and 
max 

once more the weight of the charge, 
dependent constants (Table II), and 
charge in feet . 




( eqn 2.7) 

t in seconds. W is 
Kg and are explosive 
D is the depth of the 



Table II 

Explosive Constants II 



Explosive HBX-1 TNT PENTOLITE 

Kg 14.14 12.67 12.88 

K 4 4.761 4.268 4.339 

A second and distinct phenomenon is bubble migration; so 
called because of the tendency of the bubble to migrate 
toward any nearby mass. When the bubble contacts the mass, 
it collapses and begins a series of rapid pulsations. The 
net effect is similar to that of the already described 
bubble pulse; a series of reloadings upon the structure. 

The rapidity of the pulsations of the collapsed bubble and 
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their increased localization tend, however, to magnify their 
severity . 

D. APPLICATIONS OF THEORY TO EXPERIMENT DESIGN 

Developing a basic understanding of the primary 
phenomena involved in underwater explosion theory is 
crucial to developing a successful underwater test program 
for many reasons. Two are of particular interest here, 
however : 

First, and most obviously, equations 2.1 through 2.3 
allow initial estimates of the peak pressure and pressure 
history to which a test structure is to be exposed. These 
are then input to the computer code so that preshot approxi- 
mations of the structural responses may be made; thereby 
creating guidelines for calibrating gages, determining 
charge weights, etc. 

Second, the design of the test itself is a direct conse- 
quence of the phenomena being observed. As is readily 
apparent from the earlier sections of this chapter, the 
loadings imposed upon a structure by an underwater explosion 
can be extremely complex. The incident wave may be a 
simple exponential decay dependent only upon charge and 
standoff, but the many secondary effects are subject to 
variables which cannot be readily controlled experimentally. 

The contributions of the secondary effects are signifi- 
cant and must eventually be studied. However, during these 
early phases of investigation into underwater explosion 
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effects a determined effort has been made to simplify the 
problem as much as possible. 

It has therefore been required that : 

1. The model to be studied be as basic as is consis- 
tent with stated objectives. 

2. Study be limited only to the incident shock wave. 

These two requirements are interrelated and are tied 

intrinsically to the nature of underwater explosion 
phenomena. As has been seen earlier, the incident shock may 
be approximated as a plane wave emanating from the charge 
and striking the panel. In the absence of secondary 
effects, the responses being studied here may thus be 
limited to those of the test panel alone. If, on the other 
hand, secondary effects are included, more than a single 
point source is involved and the interactions of any support 
or backing structures with the fluid medium must also be 
considered . 

Fortunately, for the purposes of this study, surface 
cutoff and bottom reflection are later time responses which 
can be neglected if study is confined to the first few 
milliseconds after detonation. Bubble and cavitation 
effects can be reduced greatly or eliminated through careful 
attention to the' test geometry. 

If an explosion occurs sufficiently near the surface 
that the initial bubble radius is greater than the charge 
depth, the hot gases forming the bubble will vent into the 
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atmosphere and be dissipated. This ven 
and, if occurring at the proper depths, 
bubble pulsation or migration. 

Bulk cavitation, too, is dependent 
etry. The region of cavitation tends t 
away from the line drawn between source 
(Fig. 2.6). The actual shape and dimen 
tion region will vary with charge size, 
the area at depths below the charge can 
expected to remain free of cavitation e 



ting is instantaneous 
prevents any further 

upon the test geom- 
o expand horizontally 
and surface images 
sions of the cavita- 
type, and depth, but 
generally be 
f f ects . 




Figure 2.6 Region of Bulk Cavitation 

(From Weidlinger Associates CUE Code.) 
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Bulk cavitation and bubble effects can Therefore 
reduced significantly or eliminated entirely if the 
geometry is such that the charge is located directly 
the plate at a depth to cause bubble venting. 

In designing the test platform [Ref. 3: pp . 19- 

series of iterations were run using the EPSA code to 
mine a charge weight and standoff distance that woul 
produce maximum test panel deflections of approximat 

9 

plate thicknesses." From these runs it was determin 
the ideal charge weight and standoff were 8 pounds o 
9 feet . 

The maximum bubble radius is expressed as a func 
charge weight and depth (Eqn. 2.6). Since it has be 
that bubble venting occurs when the ratio of charge 

3 

maximum radius is less than 0.75, choosing a charge 
is a fairly straightforward procedure once a determi 
of the desired ratio has been made. 



be 

t est 
above 

28] a 
det er- 
d 

ely four 
ed that 
f TNT at 

tion of 
en found 
depth to 
depth 
nat ion 



2 

This value was recommended by Weidlinger Associates as 
one which experience indicated would ensure stiffener 
tripping . 

3 

Any ratio below 0.75 will cause venting. Venting by 
its very nature will duct a portion of the .explosion energy 
into the atmosphere, thereby decreasing the maximum pressure 
of the shock wave. As a rule of thumb, this decrease in 
usually considered negligible until the ratio drops below 
0.50. 



34 



By selecting a depth to radius ratio of 0.50 it was 
readily determined that a charge depth of four feet was 
appropriate. Applying this to Eqn . 2.6 yields 



which assures venting. 

From the foregoing it can be seen that a test configura- 
tion that will have the desired consequences of reducing 
bubble and bulk cavitation effects can be designed. For an 
eight pound TNT charge, this platform will have the charge 
centered vertically at a standoff distance of nine feet 
above the plate. The charge itself will be located at a 
depth of four feet below that surface. The actual design of 
the test platform will be discussed more fully in Chapter V. 




7.60 ft 
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III. THE FLAT PLATE MODEL 



As stared in the introduction, the objective of this 
study has been to test the applicability of a computer code, 
USA-STAGS, to predicting the effects of an underwater explo- 
sion upon ships' hulls. To do this, it is clearly necessary 

that some basis for modelling simple sections of a typical 

4 

ship's grillage be established. 

The scaling laws used here and in the earlier work of 
this series [Ref. 3] were those developed by Dr. Raymond P. 
Daddazio of Weidlinger Associates, Inc. 

In his work [Ref. 8], Dr. Daddazio has established two 
dimensionless parameters, 3 and X, which he refers to as 
the "plate slenderness ratio," and the "longitudinal stif- 
fener slenderness ratio," respectively. These are defined 



3 



( eqn 3.1) 



t E 



and 



X 




( eqn 3.2) 



4 

The ship's hull or superst ructur al plating and associ- 
ated stiffeners. 
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where 



j = material vield strength 

y “ ° 

E = Young's modulus 

a = spacing between transverse frames 

b = spacing between longitudinal stiffeners 



t 


= plate thickness 








k 


= radius of gyrat 


ion of the 


longitudinal 


stiffener 




acting with an 


assumed ef 


fective width 


of plating 



b where , 
e 




Through sampling a variety of ship types, Dr. Daddazio 
found that typical values for ships' grillage fall into the 
range : 

1-.0 < 8 < 4.5 
and , 

0.15 < A < 0.9 

The intent of this series of studies has been to use one 
basic model but with a variety of stiffener types so that 
the responses of different stiffener geometries could be 
studied. 

The model established in Ref. 3 was a 3/16 inch thick 
test panel 18 inches long by 12 inches wide milled into the 
center of a two inch thick 6061-T6 aluminum blank measuring 
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33 inches by 27 inches. Two stiffeners with rectangular 
cross sections, 3/16 inch thick by one inch deep were mounted 
t ransversely across the test panel in such a manner as to be 
symmetric to both the transverse and longitudinal axes (Fig. 

3.1) . The entire assembly was machined as a single unit to 
avoid the inconsistencies of welding. A radius was cut at 
all corners and the whole polished to avoid stress concen- 
tration points. 

To apply Equation 3.1 and Equation 3.2 the two milled 
stiffeners were considered to be longitudinal stiffening 
members and the longer boundaries of the panel to be 
transverse frame members. Thus, the dimensions of the plate 
are, a = 12 inches, b = 6 inches, and t = 3/16 inches. Using 
the standard reference values for 6061-T6 aluminum, it was 
found that, for this original plate: [Ref. 3] 

6 = 2.02 

and 

X = 0.9 

which fall within the acceptable ranges. 

The initial plate used in this study was identical to 
that used in Ref. 3 with the exception that the stiffener 
cross section used was a tee rather than a rectangle (Fig. 

3.2) . 

Since the spacing of the stiffeners in this tee- 
stiffened plate was identical to that in the original plate, 
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Figure 3.1 



Test Panel with Transverse Rectangular Stiffeners 
[Ref. 3 : p . 13 ] 
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Figure 3.2 Test Panel with Transverse Tee Stiffeners 
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CM 



clearly the 3 values remained invariant. Given the 
changed geometry of the stiffener the effective radius of 
gyration of the plate /st if fener assembly (k) could be 
expected to change, and with it the value of X. 

Before that value X could be determined, however, the 
dimensions of the stiffener itself first had to be 
established . 

In an effort to ensure that stiffener responses were a 
function of stiffener design alone, the plastic moduli, Z , 
of the two stiffener cross sections were set as close' as 
possible to equal. For the original, rect angul ar-st i f f ened 
plate, Zp equaled 0.046875 in . 

By a process of iteration, a tee-stiffener with flange 
dimensions 23/32 inch by 1/8 inch, and web dimensions 25/32 
inch by 1/8 inch was fixed upon as giving an elastic modulus, 
Z^ = 0.046170 in ^ , which is within 1.5% of that for the 
rectangular stiffener (Fig. 3.3). 

Now, applying these dimensions to the grillage param- 
eter, X, we find that the neutral axis of the effective 
stiffener (Fig. 3.4) is: 



y = 



2 A . y . 

TV 1 



T 



= 0.2220 in 



The moment of inertia of the effective stiffener is 



= 21. + 2A . d . ‘ 
ZZ i ii 

= 0.092562 in 
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Neutral Axis : 



a 



1- 1 
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y = 



-A . y . 
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= 0.607747 in. 
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Plastic Section Modulus: 

Z = Sum of the moments about the 
" neutral axis of the cross 

sections above and below the 
axis 



Z = (b h ) (h + ^h ) 

p 112 1 



h, = 



b, = 



1/8 in. =0. 125000 in 
23 



h„ = 



32 

25 



in. = 0.718750 in. 
- 0.607747 



32 

= 0.173503 in. 



h 3 = 0.607747 in. 



b 0 = b_ = — in. = 0.125000 in. 
2 Jo 



Z = (0.718750) (0. 125) (0. 173503 + H - 0.125) 
P 

+ 4 (0.125) (0. 173503) 2 
+ Jj (0.125) (0.607747) 2 
= 0.046170 in 3 



Figure 3.3 Plastic Modulus of Tee-Stiffened Plate 
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Figure 3.4 Effective Cross Section of Tee-Stiffend Plate 
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which gives a radius of curvature: 




= [0.092562 in 4 /1.0256 in 2 ]^ 

= 0.300419 in 

Thus, for the tee-stiffened plate, 

, _ 12 40000 * 

- .300419,1 1Q x 10 6 

= 0.804144 

which is well within the parameter range required for 
modelling ships' grillage. 

The second requirement levied upon these test plates was 
the ability to demonstrate stiffener tripping. 

Stiffener tripping has been identified as the lateral- 
torsional instability of the stiffener as it becomes 
suddenly unstable and fails under load. It is identified by 
a characteristic warping and buckling of the stiffener. It 
is significant that the continued resistance of a stiffened 
structure to deformation under loading is dramatically 
reduced after the stiffeners have tripped. In earlier 
studies of impulsive loadings upon ships' grillage, tripping 
has been identified as "... a primary ductile failure mode 
for ship structure." [Ref. 9: p. 2], 

Previous observations have indicated that tripping is 
most readily initiated when an axially compressive loading 
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is applied to the top fibers of the stiffener. Recognizing 

that a flat plate will tend to deform away from the source 

of an explosion, the plates in this series of studies have 

been oriented with the stiffener sides outward into the 

5 

fluid medium. 

As stated earlier in this chapter, the original intent 

in this series of studies was to maintain an essentially 

fixed plate model, varying only the stiffener geometries. 

Unfortunately, the results of the first two explosion tests 

(on the original rect angular-st if fened plate [Ref. 3], and 

on the tee-stiffened plate here) were not satisfactory in 

view of the stated objectives. The two transversely mounted 

stiffeners were too short and the entire assembly too rigid 

to ef feet ively. demonstrate the desired stiffener tripping. 

The shock wave essentially "punched" through the panel 

material, shearing it at the fixed boundaries and at the 

6 

base of the stiffeners (Fig. 3.5). 

The results of the two initial tests indicated that a 
design change was in order. The third plate of the series 
(the second of this study) was a simplified model to permit 

5 

This is hardly the expected orientation on a surface 
vessel, but is common on many submarine pressure hulls. 

6 

Although not easy to see in Figure 3.5, not all of the 
damage to the tee-stiffened plate was done by the incident 
shock wave. Apparently, when the bottom reflection hit the 
plate and its backing structure, it was forced upward into 
its own supporting cables. Much of the damage to the stif- 
finers in particular seems to have been caused by the 
actions of that cable. 
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Figure 3 . 5 



Explosion Damage to Tee-Stiffened Plate 
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closer investigation of the shock effects themselves, 
without attempting either to retain a set plastic modulus 
or to approximate ships' grillage. 

The result of this simplification was a thickening of 

the base panel to £ inch to increase resistance to shear. 

The stiffener was restored to the original rectangular 

configuration but the thickness was increased to i inch 

7 

and the depth to It inches. Most important, in this 
attempt, a single, longitudinally oriented stiffener was 
used rather than the two transverse stiffeners in an effort 
to provide greater flexibility along the length of the 
stiffener (Fig. 3.6). 

Although the charge used to test the rectangular stif- 
fened plate was small pound of TNT), and the correspond- 
ing plate deflections were far too slight for any firm 
conclusions to be drawn, it can be seen in Figures 3.7 and 
3.8 that tripping action appears to have begun to manifest 
itself. The single, longitudinally oriented stiffener seems 
to be the superior model; and, while further testing is 
required at larger charge weights, its employment is tenta- 
tively recommended for future studies. 



7 

This decision was made somewhat arbitrarily. Dimen- 
sions were chosen, in part, to reduce the cost and simplify 
machining the plate. It is interesting to note, however, 
that, with 0 = 1.52 and \ = 1.1, this plate is not far 
outside the acceptable grillage range. 
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Figure 3.6 Test Panel with Longitudinal Rectangular Stiffener 
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Figure 3.7 Stil Tenor Del lections 
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Figure 3.8 Plate Del lections 



IV. THE USA-STAGS CODE 



A . GENERAL 

Structural Analysis of General Shells (STAGS) is the 
finite element code first developed by Lockheed Missiles and 
Space Company, Inc. during the 1960 's in response to an 
initial tasking by the Navy and the Air Force plus an inter- 
nally perceived need for a general structural analysis 
program. 

The current version of STAGS is STAGS C-l , a dynamic 
analysis code with broad applications to the aerospace 
industry which has also found application in ship design and 
engineering. A new version, STAGS 2000, is now under develop 
ment which will have an entirely changed data manager 
concept from that used in the C-l version. This change 
should make the new code both easier to use and far more 
economical of the user's time. STAGS 2000 will appear 
initially as a static analysis code, but will eventually 
exceed STAGS C-l's capabilities. 

STAGS C-l is supported by Computer Software Management 
and Information Center (COSMIC) under the sponsorship of the 
National Aeronautics and Space Administration (NASA). 

The USA-STAGS code is an extension of the basic STAGS 
routine to permit evaluation of impulsive loadings upon 
design structures in a fluid medium. As such, Underwater 
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Shock Analysis (USA) is a mulrimodule routine which is 
inserted as an intermediary and adjunct to the program pre- 
analyzer ( STAGS1 ) and analyzer ( STAGS 2 ) . The STAGS1 , STAGS2 
combination can stand alone for some types of general struc- 
tural analysis. As a result of its add-on nature, the USA 
code cannot be used in a stand alone mode, but it is not 
STAGS dependent. Versions have been created for use with 

g 

other finite element structural analysis codes. 



B . THEORY 

The characteristics which make USA-STAGS unique lie in 
its applicability to a wide variety of problems and, partic- 
ularly, in the individual techniques used to solve those 
problems. The broad, underlying theory behind those solu- 
tions is in itself common to most finite element codes. 

Briefly, the technique used by codes of this genre is to 
solve the differential equations of motion for a structure 
in order to determine its responses to its environment 
[Ref. 10: pp . 2-1 through 2-14]. By then applying the 

constitutive equations (Hooke's law generalized to include 
anisotropic materials), the stresses in the governing equa- 
tions of motion may be expressed in terms of strains. These 
strains may in turn be expressed as the first partials of 
the displacements with respect to the spatial coordinates 



O 

As of this writing, USA has been linked to NASTRAN , 
GENSAM, SPAR, and DIAL. 
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through application of the kinematic relationships. An 
equivalency can be drawn between the motions of the struc- 
ture between any two instants in time t^ and t^ and the 

extrema of the energy functional /^1(U - K)dt with respect 

0 

to all admissible variations of the strain energy U and the 
kinetic energy K which satisfy the geometric boundary condi- 
tions and continuity conditions. The equilibrium equations 
can then be derived from the first variation of the 
energy functional through partial integration with respect 
to time and the spatial variables. 

The variational approach provides a particularly effec- 
tive means for analyzing continuous systems because of the 
way in which it generates and accounts for natural boundary 
conditions [Ref. 11]. The disadvantage of the method is 
that it is limited to conservative systems. However, by 
applying the principle of virtual work STAGS is able to 
greatly extend the classes of problems to which it can be 
applied . 

In the second case; where inelastic deformations such as 
plastic strains or thermal expansions are encountered, STAGS 
treats them as "pseudo-loads." In this manner, by applying 
these pseudo-load terms to the force side of the equation, 
an initial estimate of the inelastic deformation values can 
be made and, if convergence is obtained, the nonconservative 
problem can be solved by iteration upon a series of 
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conservative problems [Ref. 10: 



of cases to which STAGS 
Figure 4.1. 

After being reduced 
the governing equations 
may be expressed in the 

M S X + D s ( X ) + B s (x) + 



pp . 2-7, 2-8]. The range 
may now be applied is shown in 

to functions of displacement only, 
of motion for a discretized system 
f orm : 

K sG> - ! 



where x is the vector of displacement components, M is 

the structural mass matrix, f is the vector of external 

forces, and is the (generally nonlinear) structural 

stiffness operator. The operators B and D include 

s s 

forces that are functions of structural deformation and 
deformation velocity, respectively. 

Equation 4.1 is the matrix equation that STAGS must 
solve to find structural displacements. For fluid/struc- 
tural interactions the generalized force vector becomes a 
much more complex entity [Ref. 12: pp . 2-1 through 2-12]. 

For submerged structures excited by an acoustic wave, 
the exciting forces are given by: 



f = -GA f ( P j + p g ) 



( eqn 4.2) 



where and P^ are the incident and scattered nodal pres- 
sure vectors acting upon the surface of the structure at the 
fluid-structural interface. is the diagonal area matrix 

associated with elements in the fluid mesh, and G is the 
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Figure 4.L Range of Problem Types to which STAGS is Appl icable 
[ Re I . 10 : p . 2-12] 



matrix which relates the structural and fluid surface 
forces. It can be seen from the development of the pressure 
characteristics of an underwater shock wave in Chapter 1, 
that for an underwater explosion the incident pressure 
vector Pj is known. The scattered pressure vector Pg on the 
other hand is dependent upon the configuration of the struc- 
tural surface and is therefore not only not immediately 
known but is also changing as the structure deforms. 

It is to solve this computational dilemma and thereby 
solve the combined equation of motion: 

M x + D S <X) + B S (X) + K S <X) = -GA f (Pg + P s > (eqn 4.3) 

that the USA code was created. 

The USA code finds P by utilizing the Doubly Asymptotic 

~ O 

Approximation (DAA) and a staggered solution procedure. 

The DAA is an approximate relation which approaches 
exactness at both the high-frequency (early-time) and low- 
frequency (late-time) limits to the passage of the shock 
front. It may be written: 

Mf Pg + pcA f Pg = pcM f Vg (eqn 4.4) 

where Vg is the vector of scattered-wave fluid-particle 
velocities normal to the structural surface at the fluid- 
structural interface. is the fluid mass matrix for the 

fluid mesh at the interface, and p and c are the density and 
sound velocity of the fluid, respectively. 
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When fluid and structure are excited by an acoustic 
wave, the incident and scattered fluid particle velocities 
(Vj and Vg) may be related to the structural response by the 
kinematic compat ability relation: 



G 



T 



X 



Yi + Ys 



( eqn 4.5) 



Solving for Vg and substituting the resultant into equation 
4.4 yields 

M f P s + pcA^Pg = pcM^(G^x - v j) (eqn 4.6) 

Equations 4.3 and 4.6 are now each expressed in terms of the 
single unknown Pg and a solution may be obtained by a stag- 
gered iteration process between the two equations for each 
time step. 



C. ORGANIZATION 

As has already been stated, USA-STAGS is a modular code. 
Functionally, it may be divided into five basic components. 

1. The structural preprocessor 

2. The fluid mass preprocessor 

3. The augmented matrix preprocessor 

4. The time integration processor 

5. The data postprocessor 

Each of these components is run individually in sequence 
utilizing the outputs from the previous modules. 
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1 . The Structural Preprocessor 

The structural preprocessor is the unit that 
assembles the structural mass and stiffness matrices, 
creates the finite element mesh, and otherwise provides an 
encoded description of. the structure's internal and ex- 
ternal degrees of freedom [Ref. 12: pp . 3-3, 3-4]. 

For USA-STAGS, the structural preprocessor is the 
STAGS1 portion of the STAGS code itself. 

Input for STAGS is made to each of ten functional 
subdivisions which correspond to the types of input 'that go 
into each. These subdivisions are: 

- Summary and Control Parameters 

- Computational Strategy Parameters 

- Data Tables 

- Geometry and Material 

- Discretization 

- Discrete Stiffeners* 

- Boundary Conditions 

- Loads* 

- Output Control 

- Finite Element Units* 

(Subdirectories designated with * are not required with all 
models . ) 

Input to each subdirectory is in the form of a 
series of data records consisting of real and integer 
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numerical data fields. This will be discussed in greater 
detail in the following section. 

2 . The Fluid Mass Preprocessor 

The fluid mass preprocessor (FLUMAS) constructs the 
fluid mass matrix for a structure submerged in a fluid. 
Inputs include [Ref. 12: pp . 3-4, 3-5] the global coordi- 

nates for the structural geometry (input from the structural 
preprocessor), plus 

- Fluid wet-surface mesh 

- Element definitions 

- Material properties of the fluid 

- Location of free surface 

- Model symmetries 

3 . The Augmented Matrix Preprocessor 

The augmented matrix preprocessor (AUGMAT) receives 
the output from the structural and fluid mass processors to 
create the matrices required for solution of the augmented 
equations 4.3 and 4.6. 

4 . The Time Integration Processor 

The functional heart of the USA-STAGS system, the 
time integrator processor (TIMINT) links back to and supple- 
ments the structural processor STAGS2 to implement the stag- 
gered solution technique; creating the solution data files 
over the time range specified. Input includes the data 
passed from AUGMAT plus 
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- Spatial geometry of charge and target 

- Initial incident pressure and pressure 

characteristics 

- Time integration information 

- Restart data 

TIMINT includes an internal postprocessing 
Generally, however, this function is left 
module . 

5 . The Data Postprocessor 

USA-STAGS Provides a selection of 
postprocessors: POSTP , POSTPR , and STAPL . 

be discussed in the following section. 

D. USING THE CODE 

As has already been stated, the STAGS 
in its concept; attempting to cover in a s 
many approaches to problem solution as can 
a single contextual framework. 

It is beyond the scope of this work to 
possible aspects of the code; indeed withi 
bilities of STAGS this approach may seem p 
limited as it is to those aspects clearly 
water shock analysis of flat plate models. 

STAGS is designed to be used for any o 
analysis [Ref. 13: p. 2-1]. 

- Linear stress analysis 

- Geometrically nonlinear elastic stress 



wave decay 



capability . 
to the following 

three possible 
Their uses will 
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- Inelastic stress analysis, geometrically linear or 
nonlinear 

- Bifurcation buckling analysis with linear or geometri- 
cally nonlinear stress state (elastic) 

- Small vibration analysis based on linear or geometri- 
cally nonlinear stress state (elastic) 

- Transient response analysis, linear or geometrically 
nonlinear, elastic or inelastic 

In this study, only two of those capabilities were used: 

(a) Small vibrations (in a stress free state), and (b) 
Geometrically nonlinear transient response analysis. Of 
these, the second alone is truly applicable to USA-STAGS. 

The small vibration analysis was used as a check on the code 
and input data. Since the actual inputs for the two anal- 
yses vary but slightly, both shall be addressed together. 

Before that can be done, however, the codes must somehow 
be brought into some semblance of conformity; for they are 
not normally provided as a single unit . 

Both codes are provided in tape form: STAGS from 4 

COSMIC, and USA directly from Lockheed. These are designed 
to be loaded onto the VAX 11/780 computer using the standard 
MOUNT MF : command . ^ 

It is suggested that, at least during the early stages 
of processing, each tape be assigned to a separate direc- 
tory, for many of the source files will be culled as the 

9 

See Appendix A. 
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executable files ( x . EXE ) 



are created . 



The source files 
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supplied with each code are specific to the code until the 
final linkage, and their sheer number can cause confusion if 
not approached methodically. 

The STAGS tape as provided contains all the necessary 
ingredients for the functioning codes except the executables 
themselves. These are left to be linked on the individual 
computer. The link commands for creation of these executa- 
bles are provided with the tape on the INFO. VAX file, which 
also contains a full listing of the tape's contents. The 
commands that are of immediate interest are LSTAGS1.COM, 
LSTAGS2.COM, LPOSTP.COM, and LSTAPL.COM. These should be 
copied into the file and executed in the normal manner 
for a program run from the terminal . ^ The result will be 
the executable files for the respective module. 

Executables are provided with the USA tape; however, 

since relinking may be required to adapt the code to local 

12 

filenames, link command files are also provided. The 

^On the VAX, the asterisk represents a "wildcard" which 
may be substituted for any filename, filetype, or version 
number (the complete listing for a file is: filename, 

filetype, version number) whenever a general category is 
desired or the listing is to be repeated. 

11 

This is done by simply typing command filename" and 
hitting the return key. 

12 

The USA tape also includes another particularly nice 
feature, the complete user's manual for each of the modules. 
This ensures that the most current version is available; 
and, more important, the version of the manual that is 
pertinent to that tape. A vital accompaniment since the 
tapes are continually being updated. 
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command files to be used in cpnstructing the USA executables 



are LINKFLU.COM, LINKAUG.COM, LINKUSAS.COM, and LINKPOS.COM. 
It is the LINKUSAS.COM file which actually performs the 
marriage between the codes; allowing the integration/ 
structural analysis functions of STAGS to be coupled with 
the USA fluid processor. 

Once all of the executable files have been created and 
tested in a full run of the combined code; all files but 
the executables may be deleted, thus realizing a considerabl 
increase m computational space. 

The basic strategy that will be followed here shall be 
to address USA-STAGS as a single code broken into five 
different modules (see preceding section); each of which can 
be run independently in stepping stone fashion once neces- 
sary input from earlier modules has been obtained. 

The first module to be run is STAGS1 . The actual input 
for this includes all information needed for both STAGS1 and 
STAGS2 execution in a single entry file. 

The first five data records: A-l , B-l, B-2 , B-3, and 
C— 1 are common to all STAGS runs. The inputs as set forth 
in the COSMIC STAGS users guide [Ref. 13: pp. 3-1 through 

3-11] are largely self-explanatory but a few points should 
be noted. 



13 

Total space required for the executables alone is less 
than 4000 blocks; compared to the almost 13000 blocks for 
the executables plus source codes. 
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First, IPOSTO , IPOST1 , and IPOST2 are applicable only to 
the STAGS postprocessors POSTP , and STAPL . If the user is 
interested only in the nodal displacements resulting from 
underwater shock these may be derived directly from either 
o-f the USA modules TIMINT or POSTPR . Thus IPOSTO, IPOST1 , 
and IP0ST2 may be set to zero thereby suppressing creation 
of model and solution data files ( * . MOD and FOR022.DAT) and 
realizing some savings in computer space. It is suggested 
that this be done in any case during early model development 
and while testing code inputs (as with the small vibrations 
analysis). If, however, strains are required; these are 
computed only within the STAGS postprocessors and IPOSTO 
must be set to 1. Considerable amounts of computer time and 
space can be used in processing a large model; time which 
would be wasted if this is not done. Similarly, IP0ST1 and 
IP0ST2 should be set to a non-zero value if displacement and 

stress plots are desired at the IPOSTlst and IP0ST2nd load 

. 14 

steps . 

Second, it should also be noted well that IFLU must be 
set to 1 if underwater shock analysis using USA is to be 
undertaken. This will cause the FOR003.DAT and FOR004.DAT 
files required for USA to be created. These files are not 



14 

This is required only when using STAPL. POSTP is 
capable of calculating displacements and stresses 
internally . 
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unduly lengthy, but for the sake of economy should be 
suppressed for non-USA runs. 



The C— 1 data record is meaningless for small vibration 

analysis and transient response analysis but must be 

1 

included. A value of 0. should be entered. 

The only true variation in data record sequence between 
small vibration analysis and transient response analysis 
comes in the next group of cards. The proper sequencing for 
small vibrations is D-2 , D-3, F-l ; while for transient anal- 
ysis it is D-l , E-l , E-2, F-l. These are the Strategy 
Parameters and the difference lies in the need for integra- 
tion over time in the transient model. The D-l record sets 
tolerances, limits the amount of CPU time, allows restarts, 
etc.-; while the E records control the actual time integra- 
t ion itself . ^ 
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At the time that this study was initiated it was 
understood that implicit integration using K. C. Park's 
formula ( IMPL = 0, METHOD = 4 on the E-2 record) was the 
only form of integration that should be used with USA. As 
of this writing, however, Lockheed informs that USA may now 
also be used with explicit integration. It has been 
suggested that this might provide a faster and more econom- 
ical solution. Due to time limitations that has not been 
tested here. 
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The D-2 and D-3 records are specific to modal analysis 
and deal with the number and range of eigenvalues to be 
calculated . 

Data records F through P pertain to the model itself. 

The records F through L provide information general to the 
entire model and are listed in normal alphabetical sequence 
once through only, although specific cards may be repeated 
within the sequence. M through P, plus Q and R are specific 
to the individual shell and beam components and are repeated 
in sequence onqe for each component (see examples in the 
appendix ) . 

The F-l and F-2 records provide grid point and stiffener 
summaries, respectively. The information in F-l controls 
generation of the finite element mesh. Inputs are simply 
the number of nodal rows and columns for each shell. Actual 
choice of which shell direction shall consist of rows and 
which columns is unimportant as long as that choice remains 
consistent for the local X and Y coordinates across the 
entire model (rows are lines of consistent X, columns are 
lines of consistent Y). Entry is made sequentially; number 
of rows ( I ) , number of columns(I) for each shell in order of 
input. The mesh created forms the master elements which 
will then be fitted into the structural geometry by the 
information contained in the M records. 

The F-2 record is a summary of shell stiffeners. Use of 
this record could be extremely useful in determining 
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generalized responses on large scale models, but was incon- 
sistent with the purposes of this study because of the 
inherent assumption that stiffeners do not deform. 

Record G-l establishes the compatibility conditions for 
the various shells. The instructions listed in the 
COSMIC-STAGS user's manual are straightforward, but refer- 
ence to Appendix B may prove helpful for establishing proper 
linkages. It will be noted that in each case, for a flat 
plate, side 1 is that nearest and parallel to the global Y 
axis while side 4 has a similar correspondence to the global 
X axis . 

Records I, J, K, and L constitute the Data Table section 
of STAGS. Again the user's manual is self-explanatory and 
will not be explicated except to note that the Material 
Record, X— 1 ; the Material Property Record, 1-2; and the 
Plasticity Record 1-3 (if used); are listed in sequence 
once for each of the materials enumerated in the B-3 record. 

If plastic behavior is expected in the model and the 
Whit e-Besseling theory of plasticity is chosen ( IPLST 
options 0 or 1 on the 1-1 record) especial note should be 
made of the user's manual's warnings concerning the slope of 
the line connecting points submitted for the discretized 
stress-strain curve in the 1-3 record. The change in slope 
at each point must be negative . Another point not mentioned 
in the user's manual, but which should definitely be noted, 
is that the first stress-strain data entry on the 1-3 record 
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must be consistent with the elastic modulus entered on the 



1-2 record. That is, the slope of a line drawn from the 
origin to that first data entry should be the elastic 
modulus. Also, while specifying plastic behavior 
may actually decrease the run time for the TIMINT module by 
reducing the number of iterations required for each time 
step, it can increase the required computer space by more 
than 100%. 17 

If plasticity is to be included in the model, there are 
also some changes that should be made in the K records. 

STAGS determines the plastic responses of the shell by 
integrating through the thickness. The input NLIP in the 
K-2 record establishes the number of integration points 
through each material layer. Since Simpson's rule is used, 
this number must be odd . NLIP may lie between 0 and 9 



(inclusive), but some care must be exercised in the choice, 
since the larger numbers will add to the length of the 



17 

A normal run for a model with about 300 nodes seems to 
be on the order of about two to two-and-a-hal f hours of CPU 
time for strictly elastic behavior. If the elastic range is 
exceeded but plastic behavior is not specified, these runs 
can easily continue for 15 to 20 hours as the code attempts 
repeated iterations for each time step (the user can control 
this somewhat by specifying the maximum allowable time for 
any iteration), and the full number of specified time steps 
may not be attainable. Inclusion of plasticity will again 
reduce the run time to roughly three to four hours. Output 
of stresses beyond the yieM point of the material, however, 
requires much larger solution data (POR022.DAT) files for 

both the TIMINT and POSTP runs. Typically, these have been 
on the order of 45,000 and 35,000 blocks respectively for 
the pressure levels used here. 
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TIMINT run; which, as has already been mentioned, is not 
insignificant for the plastic case. The default value for 
plasticity is 3. This may be acceptable for models where 
the membrane stress governs, but will not be so if there is 
marked bending stress. 

Another change that should be made for runs including 

plastic responses is to enter a value of 1 or 2 for the LSO 

entry of the K-2 record. If this is not done the USA-STAGS 

code will successfully complete an entire run, but the post- 

1 8 

processor, POSTP, will not print stress or strain data. 

The sequence of records M through R defines the proper- 
ties and loads on the shell units. As is noted most explic- 
itly in the user's manual, these records are read 
sequentially for each unit. All data in the sequence will 
be read before the M-l record for the next unit is read. 

The individual shell geometries and the defining parame- 
ters by which they are to be connected to other shells (that 
is the means by which they are to be oriented to the global 
geometry) are specified in record M-l. IGLOBE option 0 is 
useful if the local coordinates of the unit already corre- 
spond precisely to the global coordinates. Option 1 is 

1 8 

Actually, this may not be entirely true. In order to 
get any useful response whatsoever out of POSTP, it has been 
necessary to "trick" the code into producing the desired 
stresses and strains by running the module with an 
FOR002.DAT file created by a STAGS1 run for the elast ic 
case. This has been true whether or not LSO has been 
flagged, with no noticable difference in the results. 
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similarly intended to be a labor saving device. A caveat 
should be observed, however, that option 1 should only be 
used if the side connecting two shells is uniquely defined 
(i.e., no straight lines). It is strongly recommended that 
options 2 or 3 , which allow more user control by specifying 
boundaries through specific global reference points, be used 
in preference to option 1. 

An option not listed in the user's manual, but which 
could conceivably be useful for construction of cylindrical 
models is option 4. This allows the user to specify the 
location of a unit origin in the global coordinate system 
and then define the boundary in terms of one translation and 
one rotation about the reference point. 

It should also be noted that STAGS will permit introduc- 
tion of shell geometries not included in the eleven which 
are indigenous to the program. This is done through a user 
written subroutine which may be specified by choosing ISHELL 
option 1 on the M-l record. 

If the ISHELL = 1 option is chosen, the shell geometric 
properties normally input in record M-2 will be input in the 
user-written subroutine LAME as stated above. Otherwise, 
shell geometries are chosen from those provided by the code. 
Property values for these geometries should be input as 
specified in the figures of Appendix B. If plastic deforma- 
tions are expected, particular attention must be paid to 
inputs to the M-5 record. First, since Gaussian quadrature 
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is used to integrate across the mesh, STAGS will give erro- 
neous results (or none at all) if that mesh is not located 
on the midplane of the shell thickness. ECZ must corre- 
spondingly be set to 0, and IPLAST set equal to 1. No 
further inputs need be made on this record unless one of the 
shell units is located on a plane of symmetry. If that is 
the case, STAGS will not compute lateral deformations 
without further information about how the shell is to 
respond. This information is provided in the form of minute 
randomly generated displacements or "random imperfections" 
of the shell nodes from the midplane. To trigger input of 
these imperfections, IRAMP should be set equal to 1, and the 
appropriate entries made on the M-6 and M-7A records. 

The user's manual's instructions and flow patterns for 
records N through R are clear and it is believed that, with 
possible references to the input files in Appendix C, the user 
should have no difficulty completing the STAGS input deck. 

As stated in the section on organization of the code, 

USA is divided into four modules. Each of these may be run 
independently once the needed data files have been produced 
by the previous module. Except for the required access to 
those generated data files, each module is run entirely in 
stand-alone fashion. Since modular input information is 
brief and error statements are generally good, this permits 
debugging without affecting earlier inputs and generates a 
high degree of confidence in the modular output. 
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It has already been mentioned that setting the I FLU 
switch equal to 1 on the B-l record of the STAGS1 run will 
generate two data files FOR003.DAT and FOR004.DAT which are 
specific to USA. These are produced in addition to the 
FOR002.DAT and output files which are normally generated by 
a STAGS1 run. 

Each of the three data files FOR002.DAT, FOR003.DAT, 
and FOR004.DAT must be renamed after the STAGS1 run for use 
in USA. The first, FOR002.DAT, should simply be renamed for 
the file subdirectory containing the USA files (e.g., 

( * . USA )FOR002 . DAT ) . FOR003.DAT and FOR004.DAT must be 
renamed to correspond to the GRDNAM and MASNAM files speci- 
fied in the FLUMAS and AUGMAT modules, respectively. This 
may be done by entering a simple RENAME command on the VAX. 
The user will then be queried regarding the name from which 

and to which they are being changed. 

The GRDNAM file (the renamed FOR003.DAT file) contains 
the global coordinates of the structural grid points gener- 
ated for the model by STAGS1 . The FLUMAS module then uses 
this information to construct a fluid mass matrix (FLUNAM) 
as well as generate the fluid mesh data (GEONAM) and trans- 
formation coordinates that relate the structural and fluid 
degrees of freedom on the wetted surface. The FLUNAM and 
GEONAM data files will then be used in the AUGMAT module to 
create a final augmented mass matrix. 
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The module AUGMAT is the component of USA which assem- 
bles the data from the FLUMAS module and from STAGS1 into 
the specific constants and arrays required for the staggered 
solution process. The input deck for this module is very 
short and experience indicates that problems with the module 
will generally be caused by input errors to either STAGS1 or 
FLUMAS . 

USA's time integrator module, TIMINT, links the results 
of the AUGMAT module which are stored in the data file 
PRENAM back to STAGS2 and then unites with STAGS2 to conduct 
a step-by-step numerical time integration of the governing 
equations for submerged structures exposed to shock waves. 
The time integration information input to TIMINT will over- 
write earlier entries to the STAGS deck. It is suggested, 
however, that these entries should be made consistent. 

Input to TIMINT includes not only time integration 
information but also the spatial coordinates for the charge 
and for the point on the structure closest to the charge 
(the standoff). These distances are made in reference to 
the origin of the global structural coordinate as defined in 
the STAGS input. One point of caution: The surface-cutoff 

provision in FLUMAS assumes a horizontal separation between 
charge and target. As the charge approaches a point verti- 
cally above the target , some of the computed values will go 
to infinity, causing the module to fail. The surface-cutoff 
function does not automatically turn itself off and the 
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FLUMAS input must be changed accordingly if a vertical 
aspect is required. 

Provision has also been made in TIMINT to print nodal 
displacements without recourse to a post -processor . It is 
recommended that this option not be exercised. A run for 
any of the other USA-STAGS modules can be made in a matter 
of minutes. Running TIMINT can easily consume several hours 
of CPU time for a moderately sized model. It is therefore 
not time or cost effective to run TIMINT each time a new set 
of displacements is needed. 

The solution data file, POSNAM, output by TIMINT is 

specific to the USA post-processor, POSTPR. If nodal 

displacements only are desired this is by far the easiest 

and most economical route to follow. It is not, however, 

the only option. USA-STAGS provides a selection of three 

possible post-processors. As mentioned, POSTPR is internal 

to USA and may be used to provide either data listings or 

graphical displays of nodal displacements. The other two, 

POST? and STAPL , are a part of the STAGS package. A second 

19 

solution data file in STAGS format, F0R022.DAT, is created 



19 

It is actually created by STAGS2 during the TIMINT 

run . 
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and may be renamed to the STAGS subdirectory for use in 

20 

either of these. 

Both POSTP and STAPL have been designed to provide not 

only displacements but also stresses and strains; POSTP in 

tabular form and STAPL in graphical on a CALCOMP plotter. 

At the time of this writing, the Naval Postgraduate School 

did not have a graphical capability on VAX/VMS so neither 

the POSTPR-graph ical nor STAPL options were run at the 

21 

Naval Postgraduate School. 

In accordance with common usage on the VAX, the sequence 
of commands required to execute any module have been entered 
into command files (*.C0M). Examples are given in Appendix 

A. 

Once the appropriate file names have been entered into 
the command file for the module to be run, the command file 
may be executed by entering command filename " for runs 
from the terminal; or "SUBMIT filename " for batch runs. A 
file definition statement (SET DEF [*.*]) is required for 
batch runs but is not necessary for runs from the terminal . 

20 

It should be noted that the TIMINT entry, NSAVER, is 
the frequency with which responses are saved for rhe USA 
postprocessor only. The frequency with which data is saved 
for the POSTP and STAPL postprocessors is actually tied to 
the restart data in NRESET. If stress and strain informa- 
tion is required at each timestep, NRESET should therefore 
be set to 1. 

2 1 

POSTPR and STAPL plots shown in later sections of this 
work were done by Lockheed. 
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E. DISCRETIZING THE MODEL 

As its name implies, the primary building unit in the 
STAGS code is the shell element. These elements may be 
constructed into shapes which are warped through almost any 
conceivable configuration in three-space and beam elements 
may be added, but the basic unit remains two-dimensional. 

The flat plate models constructed for this study, then, 
were ideally suited for discretization using such elements. 

Each plate with its associated stiffeners was capable of 
being divided into several flat shell units composed of 
square or nearly square elements. Additionally, the even 
spacing of the stiffeners allowed advantage to be taken of 
the plate symmetries to reduce the overall size of the model 
and corresponding length of the computer run. 

The transverse tee-stiffened plate established symmetry 
in two directions permitting a quarter model to be used. 

That quarter model was then divided into five shells; two 
for the base plate, and three for the stiffener. These 
shells were numbered beginning with the two shells 
comprising the base upward through the stiffener (Fig. 

4.2). This was done to permit changes to be made to the 
code input by simply adding or deleting shells. As long as 
the basic plate dimensions and symmetries remained the same, 
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Figure 4.2 Shell Units lor Tee-Si i l ienee! Ulale 



extensive rewriting could be avoided throughout the 

22 

USA-STAGS modules. 

The STAGS code will generate the finite element mesh. 

The user need only specify the number of nodal rows and 

columns. In the case of the tee stiffened plate, an element 

0.6 inches by 0.6 inches was chosen in order to keep the 

mesh as fine as possible and yet be able to place the active 

filaments of the strain gages used during the underwater 

23 

explosion testing entirely within the element. This corre- 
sponded to 11 nodal rows and 6 nodal columns in shell 
unit 1, 11 rows and 11 columns in shell unit 2, and 11 
rows and 3 columns in each of the three stiffener shells 
( Figure 4.3). 

Since the stiffener used on the second test plate was 
rectangular, three shell units sufficed to model it (Fig. 
4.4). At the time that this plate stiffener was discre- 
tized, the limit set upon the number fluid mesh elements by 

24 



the FLUMAS module was 171. 



The increased surface area 
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Unfortunately, it was not possible to carry through 
with this plan since the plate with the longitudinal rectan- 
gular stiffener completely changed material dimensions and 
permitted only one axis of symmetry. 

2 3 

STAGS computes stresses and strains at the geometric 
centroids of the elements. Strain gages should therefore be 
centered on the element. 

24 

It turns out that this is relatively easy to change, 
but it involves technical considerations that are beyond the 
scope of this work. Questions about such modification 
should be directed to Lockheed. The limit on most accounts 
has recently been increased to 400. 
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Figure 4.3 Quarter Model of Tee-Stiffened Plate 
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Figure 4.4 Shell Units for Rectangular-Sti f fened Plate 



necessitated by the use of a half-plate model instead of a 
quarter-plate set the total number of fluid elements beyond 
that number when the 0.6 inch element was used. Rectangular 
elements 1.2000 inches by 1.2857 inches were chosen for 
the base shells in order to allow strain gage placement at 
the center of each shell. Elements for the stiffener shell 
were 1.2857 inches by 0.6250 inch. This corresponded to 6 
nodal rows and 8 nodal columns for shell units 1 and 2 , 



and 3 rows and 8 columns for shell unit 3 (Fig. 4.5). 




placed at the mid-planes of the shell units. This was not 
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originally done. In an attempt to make the USA-STAGS 
output correspond directly to the strain gage outputs from 
underwater testing, the plates were modelled with the meshes 
placed on the surfaces of the plates which corresponded to 
the surface of strain gage mounting. The half-thickness of 
the plate was then input as the ECZ value on the M-5 record, 
thereby defining the true mid-plane for the code. This 
procedure worked extremely well for elastic runs. As stated 
in the previous section, however, the Gaussian quadrature 



that STAGS uses to calculate 

the mesh surface is incapable 

25 

at any but the mid-plane. 

For each of the test plat 
taken of geometric symmetries 
the finite element model. Fo 



the plastic responses across 
of handling the integrations 

es , maximum advantage has been 
in order to reduce the size of 
r the tee-stiffened plate this 
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This important point is reiterated here because it may 
not be immediately apparent from the STAGS run. The exper- 
ence gained here seems to indicate that when a negative 
eccentricity is input to the code, perfectly reasonable (but 
wrong!) strain approximations are obtained. The first indi- 
cation of error in this method for correlating STAGS 
and strain gage data arose when a run with a positive 
eccentricity was attempted. Under those conditions, STAGS 
produced no strain data at all. This inability of STAGS is 
not incapacitating, but it is inconvenient for it means that 
the data generated in the computer run will have to be 
converted to correspond to the strain gage information. In 
theory this is simply a matter of multiplying the plate 
half-thickness by the curvature calculated in the run and 
then adding or subtracting that from the mid-plane strain. 
Given the massive size of the STAGS output, however, this 
means that either special processors will have to be 
employed or a conversion program will have to be written. 
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presents no further complications since the resulting 
quarter model is non-symmet r ic . For the plate with the 
longitudinally oriented rectangular stiffener, however, the 
result of half plate modelling is that the stiffener itself 
lies on a plane of symmetry. This presents a problem for 
STAGS. Although deformations along the symmetric plane will 
be calculated, no lateral displacements will be recognized 
since, as far as STAGS is concerned, forces on either side 
are equal. This both affects the validity of the displace- 
ments that are calculated and causes a failure to indicate 
stiffener tripping. STAGS solves this problem by allowing 
the user to either input initial geometric imperfections in 
the shell reference surface or specify that the code 
generate random imperfections. As stated in the previous 
section, this is done by flagging the IWIMP and I RAMP vari- 
ables on the M-5 record and then following the flow direc- 
tions listed in the manual through the M-7A record. For the 
rectangularly stiffened plate used here, imperfections in 

p 0 

the first three harmonics of 0.001 inch were specified." 

Finally, although not strictly related to plate discre- 
tization, there is one other point not mentioned in the 
STAGS user's manual, but which is essential if the computer 
code is to successfully compute plastic deformations. If 



2 6 

See the STAGS1 inputs for the rectangular stiffened 
plate in Appendix C. 
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plasticity is to be included, rhe variable IPLST on the 1-1 

record must be set to some value between 0 and 4 to indicate 

— 

which of four plastic strain theories should be used."' A 
non-zero positive value must also be specified for the vari- 
able NESP . This NESP value specifies the number of data 
points on a stress-strain curve for the plate material which 
will be input to the 1-3 record. There are two points that 
the user's manual fails to make clear about these stress- 
strain points . First , this curve requires stress and engi - 
neering strains beyond the yield point. Second, the first 
point of the curve must correspond closely to the modulus of 

9 g 

elasticity of the material specified in the 1-2 record." 

All the material and stress-strain characteristics for the 
6061-T6 aluminum used to manufacture the test plates here 
were taken from the Military Standardization Handbook , 
Metallic Materials and Elements for Aerospace Vehicle 
St ruct ures [Ref. 14: p. 3-184]. 

27 

The default value, 0, automatically triggers use of 
the Whit e-Besseling theory, as does the set value 1. 

2 8 

That is, the slope from the origin to the first point 
specified must be the modulus of elasticity. The White- 
Besseling theory will apply approximations to attempt a 
patch-up. If the discrepancy is large, however, continuity 
between the elastic and plastic responses will be lost and 
the resulting output will be in error. 
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V. 



TESTING THE CODE 



A. MODAL TESTING 

Determination of a structure' 
frequencies can provide a first t 
code's ability to predict the dyn 
excited system. Perhaps of even 
purposes of this study, it can be 
structural data input to the code 
proper connectivity between shell 
the input dimensions or material 
gross discrepancies in the modal 
The Naval Postgraduate School 
the few institutions to have a He 
Fourier Analyzer with modal testi 
The technique used by the H-P 
similar to the solution processes 
codes themselves for it involves 
same governing equation of motion 



s mode shapes and natural 
est of a finite element 
amic responses of an 
more importance to the 
an excellent check of the 
Failure to provide 
units or inaccuracies in 
properties will show as 
analysis . 

is fortunate to be one of 
wlet t-Packard System 5451C 
ng capability . 

5451C is in many ways 
used by the finite element 
solving basically the 
[Ref . 15] . 



Mx + Cx + Ky + f = 0 



( eqn 5.1) 



Since both x and f are functions of time, the transfer 
function, h(x), may be found by taking the Laplace trans- 
form of Equation 5.1 
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B(s)X(s) = F ( s ) 



( eqn 5.2) 



where , 



B ( s ) = Ms 2 + Cs + K 



( eqn 5.3) 



and 



H(s) = B(s) 



-1 



(eqn o . 4 ) 



If the exciting forces and displacements are known, H(s) 
may be readily found by solving the equation which results 
from substituting Equation 5.4 into Equation 5.2, namely 



For an oscillating system of order n, the transfer func- 
tion will have 2n poles occurring in complex conjugate pairs. 
Each pair of poles will cause a mode of vibration in the 
structure. Each of these poles is a complex number which 
may be expressed as 



frequency for that mode. 

Once a and go, have been determined, the resonant 
K. K 

frequency (in radians per second) and the damping ratio 
C follow directly. 



H(s)F(s) = X(s) 



( eqn 5.5) 




( eqn 5.6) 



where is the damping coefficient and oo^. the natural 
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o 

-'k 



( eqn 5.7) 




( eqn 0 . 8 ) 



The exciting forces and dis 

provided to the system may be f 

The method chosen here was to e 

29 

modally tuned impulse hammer, 
attached to the plate to record 
from which displacements were f 
int egrat ion . 

The two fundamental assumpt 
[Ref . 15 ] are that : 

1. Modal frequency and dampi 
transfer function taken f 

2. Modes of vibration can be 
an elastic structure exce 
no excitation is possible 

It should therefore be poss 

interest by simply choosing a s 

tion points on the structure an 

close enough together to minimi 

repl icat ion . 

For the two plates used in 
and tee-stiffened test models), 
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model 


086B03 


30 pcb 


mode 1 


302A07 



placements that must be 
ound in a number of ways, 
xcite the plate using a 
An accelerometer was 
the resulting accelerations, 
ound by numerical 

ions of this test procedure 

ng are constants for any 
rom the structure . 

excited from anywhere on 
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31 

were chosen on the back of each plate. Care was taken to 

32 

avoid symmetric placements of the points. ~ The accelerometer 
was fixed at a single location arbitrarily chosen but spaced 
away from the clamped boundaries. 

When any location was tapped with the modal hammer, the 
resulting forces were fed back through the hammer to the 
Fourier Analyzer; as were the accelerations from the accel- 
erometer. The Fourier Analyzer then computed the displace- 
ments and solved the Laplace transform to create the plot 
shown in Figure 5.1. Each of the peaks in this representa- 
tion corresponds to the poles of the transfer function for a 
single point of excitation. By then combining the plots for 
all points, an accurate representation of the system modes 
was obtained. It was then simply a matter of analyzing each 
of the poles in descending order of magnitude to whatever 
total number of modes is of interest to the user to deter- 
mine the dominant modes of the system. 

B. UNDERWATER SHOCK TESTING 

Underwater explosion testing for all studies conducted 
in this series to date has been done at the West Coast Shock 
Test Facility (WCSF), located at the old Hunter's Point 
Naval Shipyard near San Francisco, California. This is a 



31 

The plates were placed face downward upon a solid 
marble slab. 

32 

The actual number and spacing of the excitation points 
seems to be primarily a matter of experience. 
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Figure 5.1 Typical Modal Test Plot 
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Navy activity reporting to the Supervisor for Shipbuilding, 
Conversion, and Repair, San Francisco. As the only licensed 
facility of its type on the West Coast, they have gained a 
considerable amount of experience in conducting underwater 
explosion testing; although most of it has been large scale 
work in support of the Navy's shock qualification program 
rather than experimentation. 

1 . The Test Platform 

In order to provide the air backed condition that 
was required to simulate a ship's hull and to ensure rigidly 
clamped boundary conditions, the test plates used in this 
series of studies were securely bolted to a heavy steel 
backing structure (Fig. 5.2). This structure consisted of a 
box constructed of 1J inch thick A6 structural steel. A 
wide flange was welded around the open side of the box to 
which the edge of the plate was bolted. An O-ring -gasket 
was fit into a channel machined into the surface of the 
flange to provide watertight integrity. Access to instru- 
mentation on the inside of the box was allowed by a subma- 

33 

rine style penetrator affixed through the bottom. 

To perform the actual testing, the plate and backing 
structure were suspended by heavy steel cables below two 
large pneumatic floats in the configuration developed in 

33 

For a more complete description of this backing struc- 
ture and the rationale used in designing it, see Reference 
3, pp. 88-94. 
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Figure 5.2 



Test Plate and Backing Structure 
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Chapter 2 (Fig. 5.3). The cables were cut to length to 
assure that the nearest edges of the test panel stiffeners 
were maintained at depths as close to thirteen feet below 
the surface as possible. Once the entire apparatus had been 
towed into the bay, the explosive charge was lowered into 
position nine feet above the stiffeners in line vertically 
with the center of the plate. 

2 . Instrument at ion 

In his work, which has already been cited several 
times, LT Rentz outlined the basic theory governing the test 
of the flat plate and began developing the procedures which 
will be enlarged upon here. 

Essentially, the information sought for comparison 
with the computer simulation codes being tested has been: 

1. Fluid pressures at the plate and in the freefield. 

2. Strain in the plate. 

3. Final deflections in the plate. 

Instrumentation used to gather this data was as follows: 

a. Pressure Measurements 

Three pressure gages were used in t he tests 
conducted here. Two of these were located just as they had 
been in the previous study [Ref. 3: p. 46]. One was located 

one foot above the center of the plate and the other was 
clamped to the aluminum block bolted to the side of the 
plate. The forward edge of this second pressure gage was 
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positioned so as to be on an arc of the same radius as the 
standoff distance. 

The pressure recordings that had been obtained 
from these two gages had been so subject to interference 
from the plate structure and the rigging -of the test plat- 
form that the futility'of conducting further tests without a 
true free field pressure gage was readily obvious. Accor- 
dingly, although the first two gages were retained as 
backups, a third gage was positioned on an arc equal to the 
standoff distance but removed five feet from the centerline 
of the plate (Fig. 5.4). As will be discussed in the 
following chapter, the results obtained from this gage 
during the test shot on the tee-stiffened plate were 
excellent . 



All tests in this series utilized one-quarter 

inch tourmaline crystal pressure gages rated to 10,000 psi 

34 

The response ratio, R , is a function of the decay 
constant for the charge, 0, divided by the transit time of 
the shock wave across the gage, t^ [Ref. 16]. For the 
eight pound TNT charge chosen as the standard for these 
studies , 



9 = 0.15 msec 



34 



Defined as the ratio of the apparent peak pressure to 



the actual peak pressure. 
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Figure 5.4 



Location of Free Field Pressure Gage 
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For a quarter inch gage, 



_ 0.25 in 

60,000^- 
’ sec 

= 4.167 x 10" 6 sec 



Applying the correction factor for the oil-filled protective 
boot which surrounds the crystal, 

1.5 tp = 6.25 x 10 ® sec 

Thus , 



From the chart on page 17 of Reference 16. 

R = 0.9733 
P 

This indicates that a very good correlation can be expected 
between the actual and measured pressures , and is also the 
best indication that the correct gage size is being used, 
b. Strain Measurements 

Perhaps the single most critical aspect of 
obtaining successful data from underwater explosion testing 
of the type conducted here is the choice and application of 
the strain gages. 
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For the tee-stiffened plate, twenty SR-4 type 

3 o 

FAE-25-12-S13ET strain gages were attached as shown in 
Figure 5.5. These were the same type used in Reference 3; 
and, although rated to only 2000 microstrain, it had 
been demonstrated that their range could be extended suffi- 
ciently to provide useable data over the very earliest 
increments of the shock response. As in the earlier test, 
BLH ' s SR-4 Permabond 910 was used as the fixative, and the 
attached strain gages and their wiring was covered with 

PR-1422 A-l Base Compound, manufactured by Products Research 

36 

and Chemical Corporation of Glouster City, NJ . This base 

compound was also used on the plate with the rectangular 

stiffener where 19 MX type PAHE-250BA-350EN strain gages^~ 

(Fig. 5.6) rated to approximately 200,000 microstrain were 

used in an attempt to increase the time span over which 

information could be gathered. Also, at the recommendation 

of the distributor, BLH type EPY-150 epoxy adhesive was 

3 S 

used instead of the Bondo 910. 



3 5 

Manufactured by BLH Electronics of Waltham, MA . 

3 6 

This was done to provide and smooth compliant 
covering that would keep the gages and wiring from being 
torn off with the passage of the shock wave, and to seal the 
gages from saltwater once the rig was submerged. 

37 

Manufactured by Experimental Stress Analysis Products. 

38 

In theory, the epoxy provides a stronger bond than 
does the Bondo, which is a contact-type adhesive. 
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Fiaure 5.G Strain Gage Placement 



In practice, although some useful data was 
obtained from each test, neither was a particularly 
successful source of strain information. The number of 
gages lost far exceeded the number that performed. 

In recapitulation, a number of faults were 
found for which the following corrective measures have been 
proposed . 

1. Although the MX type strain gages greatly exceed 
the strain requirements expected in these tests , 
their use is strongly recommended. Little or no 
loss of accuracy is experienced in the lower strain 
regions through use of the high-elongation gages, 
whereas mid-range gages tend to peak right at the 
strain ranges expected. 

2. On any strain gage, the weakest points are the sold- 
ering tabs to which the gage wires are attached. 

These are brittle and will break if any bending 
stress is applied. All motion should be restricted 
to the plane of the gage. This means that the strain 
gage wires must be extremely flexible so that any- 
bending motion, either during installation or as a 
result of plate deflections, is not transmitted to 

the gage tabs. Rated strain gage wire is 
39 

recommended . 



39 

On the tee-stiffened plate, 28 gage wire made specifi- 
cally for strain gages was used. This proved to be quite 
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3. Even with flexible strain gage wire, a determined 

force along the length of the wire, as would be experi 
enced during plate deformation will break the tab, 
the wire, or the solder connection between the two. 
Although a plastic medium, the PRC Base Compound 
used to protect the gages and wires was not flexible 
enough to allow the wires to expand with the plate. 

A suggested solution to this problem is to lay the 
gage wires with some s-curves to allow elongation as 
the plate deforms. The wires should not, of course, 
be otherwise affixed. A piping of Barrier WD , a 
tacky, semi-fluid substance manufactured by BLH which 



satisfactory as far as the gages themselves were concerned, 
but presented definite problems on the human side of the 
equation. The wire was so fine that it was impossible to 
work with it . It also proved to be highly vulnerable to 
damage while the test structure was being transported or 
handled. For the following shot, the attempt was made to 
use 20 gage bell wire. This withstood handling magnifi- 
cently, but was too stiff for the soldering tabs. Annealed 
22 or 24 gage wire designed specifically for use with strain 
gages appears to be the best possible choice. 



40 

This base compound seems to be an 
medium. A couple of cautions should be 
use, however. First, while the compound 
one hour of mixing base and accelerator, 
should be dedicated to stirring to ensur 
is thoroughly mixed. Although this obse 
completely confirmed, it appears that in 
leaves a residual which can permeate wir 
actually destroy very fine copper wires, 
compound is not compatible with the Gage 
normally painted over strain gages as a 
base compound is to be applied directly 
Gage Cote should not be used. 
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Second, the 
Cote lacquer 
sealant . I f t he 
to the strain gages, 
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looks somewhat like a combination of beeswax and 



petroleum jelly, could then be applied over both gages 
41 

and wires. This piping of Barrier WD should then be 

covered by a layer of aluminum foil and then the whole 

overlayed by the usual PBC Base Compound. The net 

effect is an armored "tube" of Barrier WD through which 

the wires can run when the plate deforms. 

The fixative used to attach the strain gages should, 

itself, be a subject for careful consideration. Even 

though the epoxy is supposed to be the stronger medium, 

experience gained here indicates that it is not as 

satisfactory as the Bondo for explosion testing. The 

postulated reason for this seeming discrepancy is that 

the epoxy places a thicker layer of adhesive between 

gage and plate material. Given the extremely rapid 

loadings experienced, intense shears are created in the 

thicker substance and the epoxy literally cracks and 

spalls off. The Bondo 910 therefore seems to be the 

42 

preferable adhesive. 



41 



A bakers pastry tube works very nicely to do this. 



42 

No failures due to loading were noted in the strain 
gages attached with Bondo. Therefore, Even if the cracking 

problem had not occurred in the epoxy, the greater separation 
between gage and plate material that the epoxy creates would 
argue in favor of the Bondo. 
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Plate Deflections 



c . 

Plate deflections were by far the most reliable 
source of data that could be obtained as a result of testing 
since no electronics were involved. To get these, precise 
measurements were made at specified points on the plate both 
before and after the test shot. A jig was devised which 



took 


advantage of the 


rigid 


plate bo 


undar 


ies to 


provide 


a 


reference table 


for a 


dial 


indicator 


dept 


h gage 


( Fig . 5 . 


7) . 


All i 


measurements 


were 


made 


while the 


test 


plate 


was bolt 


ed 


onto 


the backing 


structure 


to ensure 


t hat 


plate 


warpage 


did 



43 

not affect the readings. 

C. DATE RETRIEVAL 

Signals from the gages mounted on and about the test 
plate were fed back through shielded cable to a common 
floating terminal box. From there, they were passed to an 
instrumentation shack located on one of the jetties defining 
the small bay in which the testing was conducted. A sche- 
matic of the setup is shown in Figure. 5.S. Calibration 
levels were applied to all pressure signals by Endevco 
Model 4470 Signal Conditioners. These same signal condi- 
tioners were used to provide calibration for all but twelve 

43 

Deflections have already been given for the rectan- 
gular stiffened plate in Figures 3.7 and 3.8. No post-shot 
measurements were made of the tee-stiffened plate because of 
the heavy damage which was inflicted when the plate was hit 
by one of the supporting cables. 
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Figure 5.7 Jig Used to Measure Plate Deflections 
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of the strain gage channels. Those remaining twelve 
channels were calibrated by two six-channel Honeywell .Model 
B2-6 Universal Bridge Balance units. The strain gage signal 
were then amplified by Bell and Howell CEC 1-168 DC ampli- 
fiers before final recording. All signals were recorded 
on two twelve channel Amp.ex FR 1300 tape recorders. 

Calibration levels for the tourmaline crystal pressure 
gages were all set at the 10,000 psi theoretical maximum. 

For the strain gages, however, attempts were made to use the 
computer codes to provide estimates of expected strains 
using the pressure and decay constant calculated from 
Equations 2.2 and 2.3. This was especially critical for 
the strain levels for the tee-stiffened plate where some of 
the strain gages were expected to be operating at or beyond 
their rated limits. 

Unfortunately, difficulties encountered in learning 

to use the USA-STAGS code prevented runs from being -made 

44 

before the shot dates. 

Advantage could have been taken of earlier estimates 
used for the plate tested in Reference 3 if the charge 
weights had remained the same; however, such was not the 
case . 

Since 8 pound TNT charges are listed as Navy standard 
stock items and order had been placed well in advance, no 

44 

In fact no successful runs were made until long after 
all experimentation had been completed. 
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problems were anticipated in obtaining charges of this size. 
As the date of the first test neared, however, it became 
obvious that no 8 pound charges would be available. The 
nearest alternate that could be provided in t ime for the 
scheduled test was a 94 pound TNT charge. 

It can perhaps be appreciated that dimensional changes 
to the test platform are not simply a matter of moving the 
charge location. As charge weights are changed, the stand- 
off distance and the charge depth both must be changed; 

which involves cutting new supporting cables and recali- 
brating flotation depths. All of which is extremely time 
consuming. When the personnel at WCSF were contacted about 
the possibility of changing the configuration for a 94- 
pound charge, they advised that it could not be done in time 
for the scheduled tests. 

The choice for the first test was clear, either run the 
test with the existing platform and a 94 pound charge, or 
abort the experiment. It was determined that more benefit 
could be gained from conducting the test. 

The values for P Q and 6 found from Equations 2.2 and 
2.3 for a 94 pound charge are: 

Pq = 4081.7 psi, 9 = 0.161 msec 

This initial pressure was markedly higher than that for 
an 8 pound charge. The maximum bubble radius, too, was 
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much larger. From Equation 2.6, A = 14.03 ft. This 

max 

gives a charge depth to bubble radius ratio of 



4 ft 

14.03 ft 



0.285 



which implies that some pressure release will take place 
into the atmostphere. 

For the second test shot, the problem was quite 
different. As was discussed in Chapter III, the tee- 
stiffened model appeared too rigid to successfully demon- 
strate stiffener tripping. The plate with the longitudinal 
rectangular stiffener was a largely tentative attempt to 
produce that action. It was desired to view the responses 
of the stiffener under a variety of loadings, beginning with 
small charges and concluding with ones in the 8 to 9^ pound 
range. Again use was made of what was available. In this' 
case, the available charge weight was ^ pound. 

For a i pound charge, the maximum bubble radius is 3.02 
feet; which leads to a charge depth to bubble radius ratio 
of 1.33. This implies that bubble pulsation will be a 
problem. Fortunately, however, the time of the pulse is 
relatively late, t = 0.17 seconds. Thus, by limiting the 
period of study to early times the bubble pulse can be 
effectively neglected. 

Calibration levels for the test shot on the tee- 
stiffened plate were determined by making an EPSA run on the 
plate model used in Reference 3 but with a 9A pound charge 
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weight . 



As has already been stated, 



geometric similarities 



between these two plates permitted broad generalizations to 
be made between them even if exact comparisons could not be 
drawn. The results of that run are shown in Table III. 

No similar comparisons could be made with the plate 
with the longitudinal rectangular stiffener. Strain gage 
calibrations were therefore arbitrarily set at a value at 
the higher end of the strain gage range which could be 
expected to exceed all requirements. In this case, that 
value was 150,000 microstrain. 

The results of these tests will be discussed in the 
next chapter. 

Table III 





Calibration Data for 


Tee-St i f f ened 


Plate 


Gage 


EPSA Strain* 


Gage 


EPSA Stra 


i 


6210 


11 


76800 


2 


6210 


12 


7610 


3 


17600 


13 


7610 


4 


8500 


14 


8510 


5 


21700 


15 


7610 


6 


69400 


16 


21700 


7 


69400 


17 


6210 


8 


69400 


18 


6210 


9 


100000 


19 


17600 


10 


49500 


20 


58200 


♦All 


values in microstrain. 
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RESULTS AND CONCLUSIONS 



VI . 

A. MODAL COMPARISONS 

As stated in chapter four, the primary purpose of modal 
testing is to provide a check of the finite element model. 

A properly constructed model should exhibit substantially 
the same fundamental frequency as that determined by modal 
testing and should demonstrate close coincidence with other 
frequency values across the range tested. Exact concurrence 
cannot be expected since the computer code is analyzing an 
idealized model which will not be affected by the externali- 
ties inherent in the actual test. Also, the number of 
frequencies determined by the computer is a function of the 
number of iterations than can be made over a given time 
span. Certain economies in the use of computer time must 
therefore be observed. While the code might be able to 
calculate the frequencies for a large number of modes; such 
calculation will be purchased with a heavy expenditure in 
CPU time which is not really necessary for an adequate 
check . 

Accordingly, the scheme followed here has been to limit 
the number of modal frequencies computed by the code to 
those that could be determined in two runs for each plate. 
Each of these runs can be performed entirely within the 
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5 

combined STAGS1 /STAGS2 processor without reference to the 

USA add-on. The first run simply specifies a broad range 

within which the computer can search for eigenvalues; the 

second uses one of the mid-range frequencies determined from 

the first run eigenvalues as a "seed" for further iterations 

46 

in order to produce a refined product. The actual perform 
ance of these functions is controlled by the D-2 and D-3 
records ofthe STAGS input deck. Since the procedure is 
outlined clearly in Reference 13 and example input files are 
provided for both the tee-stiffened and rectangular stiff- 
ened plates in Appendix C, further elaboration will not be 
provided here . 

Modal testing of the milled plates was conducted in the 
manner described in Section D of Chapter IV. The results of 
both that testing and the STAGS2 estimates are provided for 
comparison in Tables IV and V tee-stiffened and rectangular- 
stiffened plates, respectively. 

It will be noted that, for each of these plates, the 
STAGS2 predicted values seem to correspond very well with 

45 

Generally referred to here as a STAGS2 run. 

46 

Since, in this case, modal testing of the milled 
plates was conducted prior to the computer runs, essentially 
the same frequency range was used for the computer runs as 
for the testing. Over this range, the first STAGS2 run will 
compute five modes with their corresponding natural frequen- 
cies. Although the user's manual states that calculation of 
up to twenty eigenvalues can be specified, and that all 
existing eigenvalues up to that number will be computed, 
five seems to be the true maximum. No increase in the 
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Table IV 



Modal Comparisons for Tee-Stiffened Plate 



Measured Resonant 
Frequency (Hz-) 

629.7410 
653.2053 
721 . 1747 
772 . 3517 
. 793.0253 
1022.4708 
1310.2979 
1447.0071 
1448.4597 
1480 . 9890 
1499.5906 
1551.0845 
1556 . 0732 
1590.2947 
1747.9771 
1994.4626 
2143.0918 
2236.9248 
2259.6445 
2275.5000 
2430.0249 
2472.2305 
2571 . 7832 



STAGS2 Predicted 
Modal Frequency (Hz) 



772.7611 

957.8445 



2043.4048 

2200.9004 



2629 . 4729 
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Modal Comparisons for Rect angul ar-St if fened Plate 



Measured Resonant ’ STAGS2 Predicted 

Frequency (Hz) Modal Frequecy (Hz) 



516 . 5652 




804.2715 
924.8018 
980.5009 
981 . 1315 
1001.6130 


790 . 0o8 / 


1118.0952 
1204 . 4885 
1229.9939 
1365 . 3230 
1428 .2920 


1115 . 3547 


1651.2498 
1651 . 2747 




1852 . 1162 


1812.0698 


2040 . 0291 


2048.0015 


2154 . 0908 
2368 . 5132 
2603 . 9736 
2725 . 8828 


2176.5366 



3039.6057 



values from the modal testing, the greatest relative error 
(using the modal test value as a base) being 6.32" for the 
tee-stiffened plate and 2.16% for the rect angul ar-st i f f ened 
plate, and the average relative errors being 3. 16" and 1 . 12 % 
for the tee- and rect angular-st if fened plates, respectively. 

Less promising is the seeming failure of the code to 
predict the fundamental frequency for either of the plates. 
It is suggested, however, that this is not due to any errors 
in -the code or in the inputs to the code, but rather is 
attributable to the fact that the boundary conditions 
actually tested were not the cl amped-clamped conditions 
specified in the STAGS inputs. The H-P 5451C Fourier 
Analyzer appears to be quite sensitive to relatively 
small changes in the boundary conditions. When the tests 
were conducted, the plates were laid face-downward upon a 
heavy marble slab. No attempt was made to further clamp 
the heavy aluminum frame into which the test plates are 



amount of CPU time allowed for eigenvalue analysis (NSEC on 
the D-2 record) seems to affect that number. Computation of 
eigenvalues over a range without any seeding of known values 
is a slow process prone to convergence problems. It should 
therefore only be used to derive a value which can be used 
for the second run; or it can be eliminated altogether if 
one of the natural frequencies is already known. If such a 
frequency has been determined, it can be entered directly 
into the code and all frequencies closest to multiples of a 
fixed interval (SHIFT on the D-3 record) will be computed. 
This process is reasonably fast and will produce all 
existing frequencies up to the number (NEIG) specified. 
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milled, nor was it possible to ensure an absolutely flush 
mount between plate and slab. 

When the tee-stiffened plate was tested, the bolt holes 

had already been drilled in the frame. These bolt holes 

seemed to cause reflections which were picked-up by the 

Fourier Analyzer as additional modes. These may be observed 

to some degree in Table IV in the repeated tendency for a 

group of frequencies to cluster about a single value. The 

effect was far more noticeable on the oscilloscope traces 

obtained as each of the excitation points on the plate was 
47 

loaded. When all of the variables are considered, it is 
believed that the actual correspondence between the modal 
test results and the STAGS2 prediction is quite good and 
provides adequat-e confirmation of both the code's ability 
to predict the natural frequencies of a simple structure 
and of the finite element models input to the code for this 
study . 

B. COMPARISONS WITH UNDERWATER TEST RESULTS 

As has already been stated in Chapter V, all underwater 
shock testing for this series of studies was conducted at 
the West Coast Shock Test Facility. The outputs from the 



It is strongly recommended that future modal testing 
be conducted before the bolt holes are drilled in the frame. 
An attempt might also be made to weight the frame with sand- 
bags to simulate as nearly as possible the clamped boundary 
conditions input to the code. 
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pressure and strain gages used for the tests along with 
their associated calibration and time signals were recorded 
by two Ampex FR 1300 tape recorders. After each of the test 
shots, the data tapes and one of the tape recorders were 
transported to the Naval Postgraduate School where the tapes 
were played-back, filtered, and recorded on disk on the H-P 
5451C Fourier Analyzer for later analysis and plotting. 

As thorough a pretest analysis as was possible without 
actual disassembly was conducted upon the equipment at the 
West Coast Shock Test Facility before each shot. Calibra- 
tion levels for all sensing devices were also checked 
using technicians and equipment provided by the Naval 
Postgraduate School. Despite these precautions, the quality 
of data obtained from the tests was highly problematical. 

An excellent free-field pressure history was obtained 
from the test shot conducted on the tee-stiffened plate 
(Fig. 6.1); fully compensating the difficulties encountered 
in mounting the third pressure gage in a position where it 
would not be affected by interference. As can readily be 
observed, this curve displays all of the classical exponen- 
tial decay characteristics predicted in Chapter II. Since 
this is so, the decay constant will correspond approximately 
to the time it takes for the pressure to drop to ^ (about 
.36788) of its peak value. From Figure 6.1, the peak 
incident pressure was 3125 psi. Therefore, 
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Figure 6.1 Freefield Pressure for Tee-Stiffened Plate 
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(3125 psi)(. 36788) = 1149.6 psi 



which occurs at about 1.56 milliseconds. 

These values were input to the USA-STAGS code to conduct 
a post-shot analysis for comparison to test data. 

The peak pressure and decay constant obtained from the 
test results are much smaller than those predicted by 
Equations 2.2 and 2.3. In retrospect, this is not parti- 
cularly surprising. First, it will be recalled that Equations 
2.2 and 2.3 are derived from empirical results and are there- 
fore subject to wide variations. Even more pertinent, it 
was found in Chapter V that the ratio of the charge depth to 
the maximum bubble radius was much less than the .50 
threshold for pressure loss. The decreased pressure and 
more gradual decay are exactly what would be expected from a 
charge that had vented a significant portion of its energy- 
in the first moments after detonation. 

Given the good results of the first shot, a high degree 
of optimism was established for the results of the pressure 
data from the second, the A pound charge on the plate with 
the longitudinal rectangular stiffener. That optimism was 
not justified. 

Upon playback of the data for the second shot , it became 
clear that the amplifiers from all three of the pressure 
gages had become totally saturated. This was surprising 
since the pressure gages were calibrated for 10,000 psi and 
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the peak pressure for a \ pound TNT charge is far less 
than that . 

The methods by which the two tests were conducted 
suggested a probable cause for the pressure gage saturation. 

Once all wiring was connected and the gages tested, the 
test platform was lowered into the waters of the bay. Heavy 
halyards were attached to either end of the floats, one to 
be used as an outhaul and the other as an inhaul. Using 
these, the platform was positioned in the bay a safe 
distance from the jetties before the explosive charge was 
installed and the test conducted. 

To conduct the test with the 9? pound charge, the 
platform was centered in the bay to keep it as far away from 
the jetties as possible. There was a heavy wind blowing and 
when the cable used to carry the firing current to the 
explosive charge was taken out to the platform it was blown 
well away from the platform and the instrument cables. 

Since it was separated from the other instruments, the free- 
field pressure gage cable was farthest of all from the 
detonating cable. 

When the test shot was conducted on the rectangular- 
stiffened plate, the charge was sufficiently small that the 
platform was kept much nearer pierside to make it easier to 
handle. All of the cabling was run out to the platform 
along the inhaul . 
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The voltage used to detonate the explosive charge is 
quite large (800 volts DC) and sets up a noticeable induc- 
tance on all electrical cabling in the vicinity. In fact, 
on the strain gages, which are very low impedance devices 
and therefore relatively insusceptible to inductive noise 

this was used as a convenient, if not very accurate, zero 
. 48 

time mark. 

The pressure gages, on the other hand, are extremely 
high capacitance, high impedance gages and correspondingly 
are very much affected by induced currents. It is therefore 
believed that the close proximity of the detonation cable 
was the sole cause of the saturation observed on the pres- 
sure gages during the second test. Be that as it may, no 
pressure history was obtained for the rect angular-st if f ened 
plate, making it necessary to use the peak pressure and 
decay constant calculated from Equations 2.2 and 2.3 in 
the subsequent USA-STAGS run. 



The 


lack 


of a zero time sig 


n al 


was not 


the onl 


y defi- 


c iency 


noted 


on the strain gage 


histories . 


Three 


ot hers 


should 


also 


be mentioned by way 
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warning 
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se 



conducting similar studies. 
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This was a major deficiency in the way testing was 
conducted at WCSF. No accurate zero time was provided on 
any channel. The detonation signal was used as a best 
alternative, but since that typically covered a span of 
about one millisecond it was almost useless when only the 
first few milliseconds were being studied. 
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Before the test platform was placed in the water, all of 
the strain gage signals were balanced with their calibration 
resistors, setting their signals to zero. As soon as the 
platform was placed in the water, the change in temperature 
effected a corresponding change in the resistance of the 
strain gages. This in turn changed the balance of the 
bridge. This was an ongoing problem throughout the period 
of the tests; many of the strain gages never were brought 
into proper balance. The problem was exacerbated by the 
fact that the strain gage and its associated wiring formed 
an extremely long leg of the balance bridge which was 
subject to the environment of the bay. The balancing leg 
consisted of a calibration resistor located in the instru- 
ment shack; a very different environment from that outdoors. 

The second of these three problems is one which is 
readily correctable by the experimenter if he is aware of 
it. When the test plate and backing structure were sent to 
the West Coast Shock Test Facility, all strain gage leads 
were connected to junctions in the terminal box located on 
the side of the backing structure. Connecting these leads 
through long cable runs was one of the many tasks that was 
performed by WCSF personnel between the time that the plate 
was delivered and the date of the test shot. Their experi- 
ence was primarily with the Navy's shock test program, 
wherein the magn it udes of the strain outputs were of more 
value than the sign. Therefore, little attention was paid 
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to the polarity of the strain gage connections. Thus, the 
histories recorded can be the exact inverse of the true his- 
tory. This can be especially debilitating when comparisons 
of the shapes of strain curves is of particular importance. 

Even more debilitating, but far less obvious, was the 
last of the problems noted. 

The amplifiers on the Ampex FR 1300 tape recorders used 
to record the strain signals are calibrated for a 40% devia- 
tion from a one volt rms input signal, peak-to-peak . The 
maximum signal deflection that is theoretically possible is 
about 1.4 volts in either the positive or negative direction. 
The calibration levels input to the system are based upon 
this maximum; signals greater than 1.4 volts will saturate 
the amplifier and the playback will simply be a straight 
line at the maxima. Yet on the recordings for virtually 
all of the strain gages which survived to produce intel- 
ligible records, the output was greater than the 1.4 volt 
limit, but no concurrent indication of amplifier saturation. 
Figures 6.2 and 6.3 are offered as typical representations 
of this phenomena for each of the two plates. 

A number of possible sources for the strain gage error 
were considered: 

1. The calibration signals were improperly set. 

2. The signals were being altered either as they were 
being input to or output from the Fourier Analyzer. 

3. The playback amplifiers on the tape recorders were 
out of calibration. 
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Figure 6.2 Untranslated Strain Output for Tee-Stiffened 
Plate (Gage 11 ) . 
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Figure 6.3 Untranslated Strain Output, Rectangular- 
Stiffened Plate (Gage 15). 
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4. The calibration signal was being -input to the signal 
conditioners at the wrong location. 

The first three of these postulated sources of error were 

readily discarded. The calibration signals recorded with 

the strain gage signals were precisely those set into the 

system prior to each shot; thus ruling out the possibility 

that either the calibration levels were incorrect or that 

the playback amplifiers were out of calibration. To check 

the second point, a series of calibrated sine waves were 

input to the Fourier Analyzer using the same connections and 

equipment settings used for the strain gage data recordings. 

These calibrated signals were recorded on disk and then 

output to the plotter in exactly the same manner as the 

strain gage signals had been. In no case did the Fourier 

Analyzer alter the calibration of the signal. 

Whether the remaining suggested cause was, indeed, the 
source of error or whether the actual cause is something 
entirely different cannot be determined at this point. 
Whatever the actual cause, the effect remains the same; the 
strain magnitudes derived from the two underwater tests 
cannot be used for comparison with the results predicted by 
the USA-STAGS code. 

Although the magnitudes of the test derived strain data 
are clearly unuseable , the general shape of the strain 
curves should be essentially unaffected. In an attempt to 
resurrect the comparative nature of this study, it was 
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determined that the shape of the curves from the tests 
should be contrasted with the strain curves predicted by the 
USA-STAGS code. This attempt, too, was doomed to 
f rust rat ion . 

Figures 6.4, 6.6, and 6.8 are representative samples of 

the strain data calculated by USA-STAGS functioning as a 

post-shot analyzer of the tee-stiffened plate. Figures 6.5, 

49 

6.7, and 6.9 are the corresponding strain gage outputs. 

Similarly, Figures 6.10, 6.12, and 6.14 are USA-STAGS 
outputs for the rect angular-st if f ened plate; and Figures 
6.11, 6.13, and 6.15 are their corresponding strain gage 



histories . 

While similarities exist at some points between the 

actual strain histories and the USA-STAGS predicted results, 

50 

in general there is very little correlation. This -was 
unexpected since earlier discussions with personnel at 
Lockheed had indicated that extremely high degrees of corre- 
lation had been obtained when the code results had been 
compared to with tests conducted by the Defense Nuclear 
Agency and others. 
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Unfortunately, after all of the problems encountered 
with the strain gages, there were not nearly as many of 
these as might have been desired. 
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Once again it should be emphasized that the shapes 
alone of the curves can be compared. As already explained, 
the magnitudes of the strain gage records are not correct. 
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SHELL 3, tLEMnX l (5.1) 




Figure 6.4 Tee-Stiffened Plate, Gage 2 (USA-STAGS) 
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Figure 6.5 Tee-Stiffened Plate, Gage 2 
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SHELL L ELEMuX i (5.1) 




Figure 6.6 Tee-Stiffened Plate, Gage 8 (USA-STAGS) 




Figure 6.7 Tee-Stiffened Plate, Gage 8 
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Figure 6.8 Tee-Stiffened Plate, Gage 11 (USA-STAGS) 
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Figure 6.9 Tee-Stiffened Plate, Gage 11 
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Figure 6.10 Rect angular-St if f ened Plate Gage 1 (USA-STAGS) 




Figure 6.11 Rect angular-St if fened Plate, Gage 1 
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Figure 6.12 Rect angular-St if fened Plate, 


Gage 15 (USA-STAGS) 
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Figure 6.13 Rect angular-St if fened Plate, Gage In 



SHELL 1, ELEMENT (3.4) 




Figure 6.14 Rect angul ar-St i f fened Plate, Gage 18 (USA-STAGS) 




Figure 6.15 Rect angular-St if fened Plate, Gage IS 
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Any number of possible solutions exist for this perform- 
ance, of course; including the extremes of either USA-STAGS 
or the strain histories being totally wrong. In this case, 
however, utilizing the extended capabilities of the code 
itself may have provided the most plausible solution. 

The USA postprocessor, POSTPR , has the capability of 
calculating the fluid pressures at the fluid structural 

interface. These results may then be displayed in either 

5 1 

tabular or graphical form. When this was done for the 

cases studied here, a very interesting result was obtained. 

Figures 6.16 and 6.17 represent the predicted fluid 

pressures at the center of the tee-stiffened and rectangular 

stiffened plated, respectively. As can be seen, USA-STAGS 

predicts that intense negat i ve pressures will be generated 

within a very short time after the arrival of the shock 

wave. That is, USA-STAGS is predicting massive local 

52 

cavitation across the faces of these plates. 

Viewed in that light, the reason for the seemingly poor 
performance of USA-STAGS in predicting the strain responses 

5 1 

As noted earlier, no graphical capability exists on 
the VAX/VMS system at the Naval Postgraduate School . These 
post-shot analysis problems were run at Lockheed Research, 
Palo Alto. It was they who most generously provided the 
accompanying plots. 

52 

In actuality, local cavitation would be marked by near 
zero pressures as the water flashed into steam to fill the 
void created by the sudden displacement of the plate. Since 
USA-STAGS assumes a homogeneous medium this must necessarily 
be represented by negative pressures. 
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Figure 6.16 POSTPR Predicated Fluid Pressures for Tee-Plate 
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Figure 6.17 POSTPR Predicted Fluid Pressures for Rectangular- 

Stiffened Plate 
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can be better understood, and even the strain gage responses 
themselves become more understandable. 

For example, in Figure 6.11 an initial compressive strain 
can be observed as the plate reacts to the incident pres- 
sure . This does not occur immediately but is delayed some- 
what as the inertia of the plate is overcome. In this 
instance, the initial compressive strain is mimicked almost 
exactly by the response predicted by USA-STAGS. The strains 
then return to zero or become tensile as cavitation occurs. 

As would be expected, the responses predicted by the code 
are much more violent then those seen in the strain gage 
histories. This is so because they are activated by large 
values of assumed pressure rather than the actual zero pres- 
sure that would occur with cavitation. 

The other plots seem to confirm this basic trend: 

Initial zero or compressive strain depending upon the loca- 
tion of the element /st rain gage followed by a positive trend 

53 

in the strain history as cavitation develops. 

The actual mechanisms involved in the results that seem 
to be indicated here are not fully understood. The implica- 
tion, however, is that the velocity of the plate as it is 
being deformed is much greater than that of the fluid. A 

53 

It is suggested that Figures 6.9 and 6.15 also adhere 
to this basic pattern, but are subject to the polarity 
problem mentioned earlier. It is believed that these plots 
should actually be inverted. 
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discontinuity is therefore being created in the fluid/ 
structural interface. This discontinuity is manifested by 
the pressure drop which causes the local cavitation. 

The strain gage histories seem to imply that local cavi- 
tation is indeed occurring. Whether that is in fact the 
case, or whether the cavitation predicted by USA-STAGS is 
simply a function of the way that the plates are modelled is 
something that must still be determined, however. 

In these studies, the plate alone has been modelled. 

The structure that is "known" to the computer code can 
therefore be conceived as infinite in extent and must neces- 
sarily bear the full brunt of the pressure wave. In actu- 
ality, this is not so. The test plate is finite in size and 



pressure 


re 


lief will 


occur 


around its boundaries. 


Since the 


breadt h 


of 


the shock 


front 


is 


on the order of five 


feex for 


charges 


of 


the weight 


used 


in 


these tests (this is 
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speed of 
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decay, it will vary with charge weight), the entire assembly 

will be enveloped by the pressure front. This means that, 

almost instantaneously after the incident pressure is felt 

by the test plate, the pressure will begin to relieve around 

the edge of the plate and a negative pressure will be 

experienced on the backs of the flanges, and, soon after, on 

54 

the rear of the backing structure. Therefore the entire 



54 

Complex interference patterns can also be expected. 



137 



force of the incident pressure wave may not , in fact , be 
impinging upon the test plate; since, while that wave is 
"pushing" on the front side of the structure , it is also 
"pulling" on the back. 

Whether this decrease in force would be sufficient to 
eliminate cavitation is an open question. More testing must 
be done. For the USA-STAGS code to be able to predict this, 
however, it is not enough to model only the test panel. The 
entire assembly of test panel, milled edge, and backing 
structure must be included. If local cavitation is still 
predicted, then additional action will be required because 
Doubly Asymptotic Approximation which is the basis for 
the USA code is not valid where discontinuities such as 
cavitation exist between the fluid and the structural 
surface^ [Ref. 17]. 

C. CONCLUSIONS AND SUGGESTIONS FOR FUTURE STUDY 

Clearly, problems encountered both in conducting under- 
water testing and in using the USA-STAGS code have precluded 
any firm conclusions from being made about the predictive 
capabilities of the code. Thus, the primary objective of 
this study has not been met. It is in the fulfillment of 

55 

Lockheed has developed a processor, USA-STAGS CFA 
(Cavitated Fluid Analyzer) which can be used to model the 
cavitated region. The process that must be followed is to 
model layers of fluid at increasing intervals from the plate 
until no further cavitation has been found. At that point, 
the standard DAA approximation may be applied. 
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the second objective presented in the introduction, that of 
providing guidelines for future study, that definite contri- 
butions can be made. 

In the realm of underwater explosion testing, some 
specific problem areas were noted in the preceding chapters. 
Solutions to some of these (e.g., strain gage choice and 
attachment) were presented in the text of those chapters. 
Others of these problems are of a more fundamental nature 
and must be discussed separately. 

First, a satisfactory plan of testing cannot be devel- 
oped without an adequate source of appropriately sized 
explosive charges. The inability to match charge weight to 
test platform dimensions makes it impossible to use any 
computer code in its preshot capacity and impacts negatively 
upon the reproduceabil it y of test results. Either the 
charge weight must be standardized or the test platform 
dimensions must be altered to accommodate changes in charge 
weight. In either case, the most apparent solution seems to 
be a carefully orchestrated test plan with long lead times 
before each test shot. If a match is not possible, it is 
suggested that the shot should probably be aborted. That 
valuable information was obtained from such mismatched 
conditions in this study is considered to be more a charac- 
teristic of a program in its developmental stages than of 
well-conceived experimental practice. Until more experience 
is gained, it appears that the rationale dictating use of 
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the eight pound TNT charge as the standard [Ref. 3] is still 
valid. It is therefore recommended that that weight charge 
be maintained if at all possible. 

Second, the lack of adequate strain gage data has been 
an incapacitating flaw in the test results. Whether the 
poor quality of that data has been the result of errors in 
the equipment set-up at the West Coast Shock Test Facility, 
the interface between that equipment and the strain gages on 
the test plate, or some other cause is uncertain. What is_ 
certain is that the present arrangement of sending the test 
plate and backing structure to WCSF a few days before the 
test shot for attachment and then allowing only a few hours 
immediately before the shot for equipment testing is 
insufficient. Development of an in-house instrumentation 
capability at the Naval Postgraduate School is considered to 
be essential. This would permit the gages and electronic 
equipment to be integrated into a single system. Thorough 
testing could be conducted before the actual date of the 
test; and, equally important, responsibility for the outcome 
of the test would reside in a single organization. 

Changes in the electronic instrumentation have already 
been recommended [Ref. 3]. Enactment of those changes would 
undoubtedly improve the data responses. A few other modifi- 
cations suggested by the results of this study are consid- 
ered equally important, however. 
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1. Rather than using calibration resistors, the strain 

gage bridge should be balanced using strain gages 
identical to those on the test plate. These should 
be attached to an aluminum plate in the terminal box. 
This has two effects: First, the balancing gages 

would be subjected to the same changes in temperature 
as the plate gages. Second, the balancing leg of 

the bridge would be of the same length as the measure- 
ment leg. Combined, these would make the bridge 
far less sensitive to environmental changes. ' A 
zero balance could be maintained and calibration 
fluctuations would be less likely. 

2. A zero-time triggering signal should be superimposed 
upon all instrument channels. This would permit an 
exact assessment of the time of initial plate response 
to the incident wave. 

3. Further shielding and grounding of the pressure gage 
cables should be attempted to prevent the saturation of 
the pressure gage amplifiers through currents induced 
by charge detonation. 

From the results that have been explored in the 
preceding section of this chapter it is clear that certain 
questions need to be raised about the validity of the flat 
plate model that has been used here. These same questions 
should be asked about the ability of USA-STAGS or any other 
code to handle that model under the conditions imposed. 
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In particular the question of whether the local cavita- 
tion predicted by USA-STAGS is actually occurring or is 
simply a result of the choice of model must be addressed. 

It is recommended that this study should be redone, but 
with the backing structure included as a part of the basic 
model. POSTPR could then be run to determine whether cavi- 
tation was still occurring across the plate. If cavitation 
was still being indicated, then either a new model would 
have to be found or the cavitated region would have to be 
modelled using a code capable of handling 'it. 

Whatever the findings of continued tests, it may have to 
be recognized that increasing the scope of the models 
studied may put USA-STAGS beyond the reach of mast ers- level 
research at the Naval Postgraduate School. USA-STAGS is a 
diverse, highly capable code which has the reputation for 
high levels of accuracy. The price that is paid for such 
capability, however, is size and complexity. Given the 
limited time and computer resources available to most 
students at NPS , USA-STAGS may be beyond their capacity. If 
such research is undertaken, it is most emphatically recom- 
mended that a large block of time, weeks, perhaps, be spent 
at Lockheed Research in Palo Alto, California learning the 



code and perfecting models. 



APPENDIX A 



SOURCE DATA 



This appendix contains the information listings supplied 
with the STAGS and USA source tapes. Section A is a full 
listing of the contents of the STAGSS4 tape along with 
instructions for mounting the tape, command files for 
creating and linking executables, and sample command files 
for running the code. Section B contains similar listings 
for the USA source tape excepting that no command files are 
listed. USA places the link commands on the tape itself for 
immediate access. Command files for running the modules are 
identical in construction to those for STAGS and can be 
adapted from those in Section A. 
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A. THE STAGSS4 SOURCE TAPE 



s £r’ E ICN Ao^UET 1 R S4 



ST AG SC 1 PROGRAM i 
VAX /VMS OPERATING St'STE*' 

1 TAPE FORMAT 

9-Track 
1600 BP I 

Labeled (STAGS4) 



TAPE CONTENTS 



lie 


Name 


Type 


Description 




1 


INFO. VAX 


Text 


Table of 


concern rs 


# procedure's 


2 


STAGSU FOR 


Sourc e 


FORTRAN 


prog "ami 


for STAGS utility lib 


3 


BLICMF MAR 


Source 


MACRO program. *cr 


STAGS utility lio:an 


4 


JOBSTA. MAR 


Source 


MACRO program for 


STAGS utility l iDrar 


5 


ST AGS 1 FOR 


Sou^c e 


FORTRAN 


programs 


f o r STAGS 1 1 ibrary 


6 


MOUNTEL. FOR 


Sourc e 


FORTRAN 


programs 


^or the MOUNT element 


6 


STAGS2 FOR 


S o ur c e 


FORTRAN 


programs 


for 5 TAGS 2 library 


9 


STAPL. FOR 


Sourc e 


FORTRAN 


programs 


for S^APL library 


9 


PQSTP FOR 


Source 


FORTRAN 


p r o g - am s 


for PQSTF library. 


10 


ST AGS 1 . 0L3 


L l 0 r a r y 


L i b r a r y 


of S T a C- 5 1 


programs 


1 1 


STAGS2 0L3 


Library 


Library 


of S"aOS2 


programs. 


12 


STAGSU 0L3 


L i b r a r y 


L l b r ary 


of S "aGE 


utility routines. 


13 


POSTP. OLB 


Li b r a r y 


Library 


of POSTP 


programs 


14 


STAPL. OLB 


L 1 b r a r y 


L i b r a r y 


or STAPL 


plot programs. 


15 


NICE OLB 


Library 








16 


TP FOR 


Source 


Thurston 


Processor. 


17 


TP. OBJ 


Ob jec t 


Thurston 


Processor. 


18 


2 SYS. FOR 


Source 


2-system 


> utilities (for TP). 


19 


RES. FOR 


Sourc e 


R e d uc e d 


Equation 


Solver 


20 


RES OBJ 


Object 


Reduced 


Equation 


Solver 


21 


ZSYS OLB 


L i b r ary 


Z-sy s tem 


utilities object library 


:-36 


C ASE* . + . * 


List 


Input and output 


for 5 sample cases 



3. LOADING FROM TAPE 

Use the following control cards to copy all files to disc 

% MOUNT MT AO: STAG84 TAPE 

% COPY TAPE: *. * * C ] 



Files may also be read ind i dually from cape with the 
following control cards (e g , to load STAGS! GL3 ) . 

4 COPY TAPE* STAGS1. OLB C] 

Complete sample case output may be oo earned as hollows 
4 PP INT /HEAD CASE* *. * 
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LISP COM 



programs a n c c - 



(Compile 



i -Ok T A G d 0 "OR / NOL 1ST 
$ L I BP AF y / CREATE STaGSu OLB ST AG 3 U 
s MACRO BLIGMF MAR 

4 macro JOBSTA MAP ' 

t LIB ST AC SO B L I GMF , JO 3 STA 

i FOP STAGS 1 . "OR /NOLI 57 
$ L I 3 P AR V / C R E A TE STAGS! GLS S T a G S 1 
$ FOR MGUNTEL FOR /NOL 1 ST 
$ LIB STAGS 1 MGUNTEL 

t FOP STAGS2 FOR/NOLIST 
s LIBRARY /CREATE S i A G 32 OLB STAGS2 
s FOR MGUNTEL. FOR/NOLIST 
s LIB STAGS2 MQUNTEL 

* FOR STAPL FOR / NOL I ST 

* LiBFARY/CREATE STAPL OLB STAPL 

* PQSTP FOR/NQLIST 

5 LI BRARY/ CREATE POSTP OLE PQSTP 

s FOR TP FOR/NQLIST 
i FOR ZSYS F 0 P/NCLI 3 T 
S LIB/CREATE ZSYS ZSYS 

* FOR RES. F 0 R/NQLI 5 T 

B. LSTAGS1.COM (Load ST AGS 1 ) 

% L I NK /MAP = ST AGS 1 / BP I EF / EX ECU TABLE =ST AGS 1 - 
STAG51 QL3/LIB/ I NCLUDE=MA I N 1 , - 
STAGSU. OLB /L I B 

C. LSTAGS2.COM (Load STAGS2) 

s LINK/MAP=5TAGS2/ BRIEF /EXECUTABLE-STAGS2- 
STAGS2 QL3/LIE/ I NCLUDE=MA I N2, - 
STAGSU OLS/LIB 

D. L 5 TAPL.COM (Load STAPL) 

S LINK/MAP= 5 TAPL/BR I EF/ EXECUTABLE= 5 TAPL- 
STAPL. 0 LB/LI 5 /INCLUDE= 3 TAPL, - 
STAGSU. OLB /LIB. - 

(File containing local plot routines) 

E. LPOSTP COM (Load PQSTP ) 

% LINK/riAP=POSTP/BP IEF /EXECUT ABLE=POSTP- 
POSTP OLS/LIB/ 1 NCLUuE=P05TP , - 
ST AGS 1 OL8/LI9, STAGS2 OLS/LIB. - 
STAGSU GLB/LIB. NICE. 0L3/LIB 
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i •* 3 S I G N C A£Ei IN? F 0 h 00 5 
i ASSIGN C a £ E i GOT FGRCOc 
4 ASSIGN C ? £ £ 1 MOD FOR 021 

4 RUN £ T A G £ 1 EXE 

5 PRINT C A5E 1 OUT 

*: G. STAGS. COn ** - (Execute STAGS 1 anc STAGES) 

% ASSIGN CASE 1 INF F0PC05 
$ ASSIGN CASE!. OUT FQROGS 

4 ASSIGN CASE! MOD FOR 021 (optional if msdsl plots oesired ) 

$ A UN STAGS 1 EXE 

4 ASSIGN CASEl RS T F0RC2G (previous solution data* r eqd for restart 
4 ASSIGN CASEl SOD FuF G 22 (optional if solution data snoulc oe saved 
4 RUN ST AG 3 2 EXE 
4 PRINT CASEl GUT 

H. PQSTP COM (Execute STAGS1 and POST? ) 

4 ASSIGN CASE A INF F0RG05 
$ ASSIGN C ASEA OUT FGP0G6 
4 RUN STAGS 1 . EXE 
4 ASSIGN CA3E4 PPN F0R005 
4 ASSIGN C ASE a POP FOR 006 
4 ASSIGN CASE4 RST FGR020 
4 ASSIGN CASE4 SOD FGRG22 
4 RUN PQSTP EXE 
4 DEL FGR009 DAT 

I. STAPL. COM (Execute STAPL) 

4 ASS I GT4 CASEl PIN F0R005 
s ASSIGN CASEl RST FGF020 
4 ASSIGN CASEl MOD F0RG21 
4 ASSIGN CASEl OUT F0R006 
% RUN STAPL. EXE 

4 (Execute plot file on local plot device) 



(Mag be omitted FOR 002 DAT saved) 

(previous solution data from CASE4 SOD) 
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B. THE USA SOURCE TAPE 



. EP 5 tC-’I r U[\i a _ . E i ^E ►- 

vaa Tape load i ng ine^^^'Ctione 



' r i : ^oe c c n t a 1 n s s y m o o i i : s . me iuc- : ' o c ■? a u r e s • : i j ? c r i : o - *] r i ? * 

« r o : u t 9 c 1 9 5 . ana manuals -or a ni *0 9 r of U E -• •: r o c 9 i iO r s All f l i e i r a , ? 
D99P :r?<3t9C with COPY commands and can oe "^nevoc in trie same manner 
with only a lines of inpu: a c the - e r m i r j 1 Examination of tn? tao? 

directc”y will reveal that :n? files nave * a ^ i o u s cycle numbers ana uni? 
13 the * e y to Keeping tnose files together - a t oelong to a particular 
processor. The following t a 0 I e contains t n i s identification code 



USA PROCESSOR 



C > LLa 



FLUMAS 
AUGMAT 
TIM INI 

post=p 
MGDL 1 2 

To load the TIMINT files, for example, go r : ■ n e directory jn?r? they 

i s 

This 



ine J cnaracter extension- on tne r : » e names denotes the tsipe or 
file according to the following table 

EXTENSION TYPE 



are to reside and simply do a COPY TAPE * % '3 C] command where TAPE 
the internal name that has oeen assigned ft t S o magnetic tape tjr.it 
process is simply repeated until all aesireo directories h a * e Oeen 
loaded 







. POP 
. INC 
. VAX 
. OLB 
. COM 
. EXE 
. UAL 



INCLUDE PROCEDURES 

INCLUDE PROCEDURES 

SYMBOLICS 

OBJECT LIBRARIES 

LINK STATEMENTS 

EXECUTABLES 

MANUALS 



For example. 
LINKUSAS. COM 
that all LINK 
referenced in 
contain users 



AUGMAN UAL is a file containing 
is the LINK statement for USA-S 
statements must De modified si 
them reflect names that are m 
names that will not exist at o 



the AUGMAT manual while 
TAGS It should o? noted 
ightly as the directories 
use here at Lockheed and 
ther VAX installations 
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I 



Dl^cCtO-y wP [ Tr 3] 








A D n ** 5 S EXE- ! 


FT/ 57 


*-SE c - 1 4 A . 


00 


AM T A r ,jp , 1 


1 / 1 


1 7 — M A y — 1 p 0 - 


vO 00 


AMGE01 r DP . 1 


1/ 1 


4 - JUN- l p 8 — 


0 00 


AMGE02 r D p < 1 


1/1 


4 — JUN" i p 9 - 


oo 


AMGEG3 r DP , 1 


i / 1 


19-JUL- ive.r 


0-; 00 


AMGEQ4 P DP / 1 


1 / 1 


4- JUN- 1 93- 


r 5 GO 


FLJL I 3 OLD. 1 


*9/69 


6-3EF- i 9S-i 


00 00 


FLUMAN UAL; 1 


37/37 


— SE p - 1 p 9 “i 


00 GO 


FLORAS FORj 1 


74/ 74 


6-SEP- 1 98- 


os: . oo 


L I NX FLU. COM, 1 


1/1 


6 — SEP “198- 


./O 00 


NAMrLU rDP; 1 


1/1 


22“APR - 1951 


vo oo 


AMGEOl PDP. 2 


1 / 1 


4— JUN- 1 98 4 


uc 00 


A M G E u 3 PDP. 2 


1/1 


22-JUL-1P85 


OC . 00 


aijGLIB OLD. 2 


25/25 


6— SEP - 1 99- 


00: 00 


AUGMAN UAL, 2 


12/12 


13-MAY- 1935 


aO: 00 


AUGMAT FOR, 2 


23/23 


&-SEP-19S- 


50 00 


FIuOAT POP; 3 


1/1 


23-MAY- l p 04 


00 00 


LINK AUG COM , 2 


1/1 


6-3EP- 195- 


uO. 00 


NAMAUG PDP, 2 


1 / 1 


22-AP R- 1 981 


00 00 


FFEPROC EXE. 2 


42/ 42 


6-SEP-190- 


00. 00 


TP AN 5 F r DP ; 2 


1/1 


29-MAt -1984 


00 00 


CH'GDAT PDP, 3 


1 / 1 


30-MAR- 193 : 


00 00 


COMADD PDP; 3 


1 / 1 


30-MAP - 1 p S A 


00 00 


CSKFAC. INC, 3 


1 / i 


29-MAR-: 9Q«- 


00 1 00 


DATFIL. PDP. 3 


1 / 1 


1 5- JUN- 1 R0- 


00 00 


DATINT PDP, 3 


1/ 1 4 


30-MAR- 1984 


00 oo 


INTLIB OL3 . 3 


52/53 


6-SEP- 1 98 J- 


00: 00 


IMTMAN UAL . 3 


22/22 


3-N0V-1983 


00 00 


LINKUSA con, 3 


1/1 


6-SEP- 1 984 


00 00 


NAM I NT . PDP; 3 


1/1 


30— MAR - 1 98 4 


00 00 


T I M I NT FOR, 3 


58/58 


6-SEP- 1984 


00 00 


UNWASH EXE. 3 


70/70 


6-SEP- 1 984 


00: 00 


CQMADD. PDP; 4 


1/1 


30— MAR— 1984 


00. 00 


DATFIL. PDP; 4 


1/1 


30— MAR - 1 984 


00- 00 


LINKPOS. COM; 4 


1/1 


30— MAR — 1 98 4 


00: 00 


NAMINT. PDP. 4 


1/1 


30-MAR— 1 984 


GO: 00 


POSLIB. OLB / 4 


38/38 


30— MAR — 198 4 


00 00 


p OSMAN UAL i 4 


15/15 


27-MAR- 198- 


00. GO 


POSFRO FOR. 4 


41/41 


29-MAR- i 9S4 


00 00 


P05TPR . EXE/ 4 


61/61 


30— MAR— 1 984 


00 : GO 


CHGDAT . PDP; 5 


1/1 


5— APR— 1983 


00 00 


COMADD. PDP. 5 


1/1 


18-MAR-19S3 


00: 00 


DATFIL. PDP. 5 


1/1 


1 5- JUN- 1 984 


00. 00 


DATINT PDP; 5 


1/1 


15-JUL-1982 


OC 00 


LINKUSAS. COM; 5 


1/1 


6-SEP— 1 984 


00 00 


MODLI3 FOR. 5 


10/10 


6-SEP-l p 04 


oo- 00 


MODLI3 OLD/ 5 


9/9 


6-SEP- 1984 


00 00 


NAMINT. PDP/ 5 


1 / 1 


22-APR-1982 


GO - 00 


USAS EXE, 5 r 


223/223 


6-SEP - 1 Q 84 


00 00 


USTLI3 OLD / 5 v 


7/7 


6-SEP-1984 


00 GO 


UST3YM FOR; 5 


11/11 


6-SEP- i 98 4 


00 00 


3 IOSUP MAR. 8 


9/9 


15-AUG-1981 


00 00 


CC FOR/- 8 


1/1 


9-SEP- 1 982 


GO 00 


CGMERR. PRC, 8 


1/1 


6-JAN— 1 Q 32 


00 00 


DMGASP FOR, 3 


13/ 18 


0- JAN- 1 p 82 


00 00 


DMGa SP 03 J. 5 


12/12 


8— JAN— 1 v8 i. 


00 00 
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J. _ _ '_,v 


P C • E 


1/ 1 


c - ...at. - ; 4 £ .. 


„ _■ J 


iO‘"CjM - 


P C , ~ 


1 / 1 


o - j an— : w £ i, 


”* - 


icmd-t = 


4. “ 


it*. 


: - ..A?,- ! == r 


: : oo 


1 1 N C Er 3 


R C , S 


1/ 1 


— — JH|'. * . ^ Z — 


- 1 0 w 


I C'-i£jV 3 


P C • s 


1 / 1 


Z. W — 1 . — ^ — 


-■ . j 


N I D E C - E 


, 5 


2 ^ 5 / 2 - 5 


1 0- * JG- 1 93 - 


' ,: I 00 


OUC n TE5 


o 


1/ 1 


a! 1 -NOV— 1 9S j 


I 00 


5 f ME I G 3 


PC, 5 


1/1 


t- JAN- 1 RSI 


00' 00 


SYMBOL 3 


PC, 5 ... “ ■ 


’ 2/2 


o-v^l'J- 1 332 


o: oo 


B'rM DAT PRC, 5 


1/1 


c- VAN- 1932 


00 00 


TE5TGA5P 


Tisi S 


1/1 


6- van- i 9g r 


OC 00 


Sr I ME 


FjDP 9 


12/12 


2- AUG- 19£2 


00 00 


SK>L INE 


0L3 , 9 


9/9 


10- AUG- 1932 


* J u 


CVL daa, 


1 1 


3/3 


30-MAR- 1 99 A 


00 CO 


CVL FLU, 


1 1 


3/3 


30-MAR - 1 9S~* 


00: 00 


CVL. GEO, 


1 1 


2/2 


30-MAR- 1 98- 


00 00 


CVL PQ5, 


1 1 


161/161 


1 - APR - 1 90 A 


00 00 


CVL. PRE, 


1 1 


5/5 


30-MAR- 193 A 


OC 00 


CVL RST, 


1 1 


20/20 


1-AP=- 1 93 A 


00 00 


CVL UEC, 


1 1 


1/ 1 


1 3- JAN- 1 939 


DC 00 


CVl USD, 


1 1 


2/2 


13- JAN- 193 A 


00 00 


C t L AIJG I 


NP; 1 1 


1 / 1 


1 7-FE3- 1 93 2 


00 00 


C -.'LAUG OUT. 1 1 


/ “3 

c. / c- 


30-MAR— 193V 


00 00 


C t Lr LU I 


NP, 1 1 


1/1 


4'— ..JAN— 1 99 P 


00 00 


C'Tt_FLU GUT, 1 1 


39/39 


30-MAR - 1 93- 


00 00 


CVL I NT INP; 1 1* 


1 / 1 


1 -APR- 1 989 


CD 00 


C V>_ I NT OUT; 1 1 ^ 


eo/ao 


1-APR-l 99 a 


CO 00 


CVL°OS INP; 11 


1 / 1 


1 -APR- 1 93- 


00 00 


CVLPOS OUT, 1 1 


113/1 13 


1 -AFR- 1 909 


OC* 00 


CVLSTG I 


NP, 1 1 


1 / 1 


9— JAN— 1939 


OC . 00 


CVLSTG OUT, 1 1 


13/13 


1 - APR - ’909 


GO . 00 


F0R002 DAT; 1 1 


120/120 


1 - Ap R - 1 90-* 


00: CO 


FQP022 DAT, 1 1 


15/15 


t -APR- 1939 


00* 00 


BEAM DAA ; 1 5 


3/3 


3— JAN— 1939 


00 • GO 


BEAM GEC 


); 15 


1/1 


3 — JAN— 1939 


00 00 


BEAM r RE; 1 5 


2/2 


3- JAN— 1989 


00: 00 


BEAM. USC , 1 5 


1/1 


3- JAN- 1 989 


GO 00 


BEAM USD; 1 5 


1/1 


3- JAN— 1 98 9 


00: 00 


BEAMAUG. 


INP; 15 


1/1 


3- JAN- 1984 


00: 00 


BEAMFLU. 


INP; 15 


1/1 


24- JAN- 1 984 


OC 00 


BEAM INE 


INP, 15 


1/1 


3- JAN- 1934 


00 00 


GEAMINS. 


INP, 1 5 


1/1 


3-JAN- 1934 


00* 00 


QEAMPSE. 


INP; 15 


1/1 


3- JAN- 1 989 


OC 00 


E E AMP SS 


INP; 15 


1/1 


3-JAN- 1934 


00 00 


BEAMSTG. 


INP, 15 


1/1 


3-JAN- 1934 


00: 00 


Total ot 


• 101 files, 


1954/1954 


blocks. 
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APPENDIX B 



BASIC SHELL TYPES FOR USE WITH STAGS 



This appendix contains a detailed listing of the eleven 
basic shell types used with the STAGS code. Drawings are 
included to show numbering conventions for edges and 
corners. This information can be particularly useful for 
discretizing the shell unit and specifying compatibilities 
with .adjoining shells. (All figures from Reference 13.) 
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x.X 



PROP(l) = X, 
PROP(2) = X‘. 
PRO?' 3) = Y7 
PROP(a) z Y, 

Subs c ripe 3 re ter :o 
corner num b e r . 




QUADRILATERAL PLATE: 

= (1-XJ • (l-Y), : 2 *X(i-Y). f 3 = XY, 




x,X 



4 

x(XY) = E f. X. 

i=i i i 



1 ' lx ' a = S, f i Y i 



PRO P( 1 ) = X. 
PRC P(2 ) = Y , 
PRO P( 3 ) = X, 

pro p(4 ) = y; 
PRO P( 5 ) = :<3 
PROP(o) = Y ; 
PROFIT) = X , 
PRO P(3 ) = Y^ 



X(l-Y) 



A X X u L A R. 




y z = r = A 
x = t(X, Y) = 0 
y = z (X, Y) = X sm Y 
c = h (X, Y) = X cos Y 

PROP(l) = R a 
PROP(Z) = R b 
?ROP(3) = 9 a 
PROP(4) = 3 k 



Lac era i disp la cement. positive up 
(negative x)! 
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CYLINDER: 




| 




y *■ z = R 
x = £ (X. Y) = X - X. 
y = g (X. Y) sRshV 
s * h (X, Y) = R cos Y 

PROP(l) = X L 
?ROP(2) = X 4 
PRO P(3 ) = Y ■ =9 
PRO?(4) = Y, = 3 h 
PROP(5) = R~ 

L = X 4 - Xl . 




i \ 



CONE: 

y~ +• z“ = r" = ((X - X. ) • :ia; - R- )" 

Can s * (R b - R a ‘) / (X 4 - X t 1 
x = £(X, Y) = X - X t 

y * g (X, Y) »((X - X t ) • can 5 - R a ) . 3 :n Y 
z = h (X, Y) = ((X- X L ) • tan 3 - R a ) - cos Y 



; \ 




PROP( 1 ) = 
PROP(Z) = 
? RO P ( 3 ) = 
PROP (4) = 
PROP( 5) = 
PROP(6) = 



R 

a 



ion) 



v 
* ? 
a: 



a 

- a 




may be > R^ 
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ISHZLL = ? 



oPriz-RZ: 




x” + y~ -r z~ = R” 

X = £(X, Y) = - R cos X 
y = g (X, Y) = R sin X sis 
z = h (X, Y) = R sis X cos 



? ROP(l) = X, = , 

1 a 

?RCP<2) = X, = , b 
PRC?(3) = Y, = 9 

1 a. 

PROP (4) = Y_ = 
PROP(5) = R 



TSKHL-L r 3 



TORUS: 



x = t(X, Y) = -R, cos X 




PRO?(l) = X, * 7 



a 

J b 

3 a 

9 b 



PROP( 5) = R 
PROP(S) =■ R. 
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C0:;z/ CYLIXDZR: 





1 




x,X 



2 1 za 

tan 3 = (R , - R )/ (X , -X. ) 

zo za 41 

r = (X - X, ) • tan y = tar, c • R /R 

y l ya za 

x = f(X. Y) = X - X L 

y = g(X.Y) = - m Y e 

z = h(X. Y) = cos Y * 

Y e = arctaa (tan Y • R /R ) 

za ya 

PROP(l) = X, 

PROP(2) = XI 

? RO P ( 3 ) = Y^ = 9 

P RO P ( 4 ) = Y, = 9u 
P RO P ( 5 ) = R“ 

PROP(6) = RR* 

P RO P ( 7 ) = ?4 b 

R mav be > R , . 
za ' zb 



I5KZLL = 10 




PARABOLOID : 

7 ? I 

y * a r" = 4R * (x + R, ) 
x = £( X, Y) = X - X t 
7 = -v /R a (X R s ) • sin V 
2 “ 2,yR a (X ' R b ) * cos Y 

PRDP(l) = X. 

PRO P( 2 ) = X 4 
PROP(3) = Y t = 

P RO P ( 4 ) : Y) s 9 b 
PRO P( 5 ) = R a 
P RO P ( 6 ) a R b 
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ISHZL- = LI | 




zzzipscio . 
? *> 



Y z 




= L 



R- 




X - a reran (ran X ' P„ / R.^) 



PRC?(5) r R 

x 

PRCP(6 ) = R 

y= 

HYPERBOLOID : 

7 7 



R w 

a 



x = X-X 



(x - RJ 

_ i 

Rr 

o 



= l 



y = (R /R. ) VRT - ( X - X i - R ) 
a o □ i c 



(R /R. )V.v - (X-X, -R P - cos Y 
a o o i c 



prop(i) = :< L 

PRO P( 2 ) = X 4 



?R0P(3) = Y 



PRO P( 4 ) = Y, 
PRO P( 5 ) = R 

a 

PR OP( 6 ) = R. 

o 

PRCP(T) = R 
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APPENDIX C 



USA-STAGS INPUT FILES 

This appendix contains the input files used for the 
complete USA-STAGS runs for both the tee-stiffened and 
rect angul ar-st if f ened plates. These are listed in the order 
in which they would normally be run: 

- STAGS2 (Modal Analysis) Input File 

- STAGS1 Input File 

- Fluid Mass Input 

- Augmented Matrix Input 

- Time Integrator Input 

- USA postprocessor (POSTPR) Input 

The input files for the tee-stiffened plate include an 
example of STAPL input. The files for the rectangular- 
stiffened plate contain an example of POSTP input file. 

Since these last files are readily generalized, it was not 
considered necessary to provide examples for each case. 
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TEESTAG32.INP; 10 ? - 3 E 3 - 1 ? 5 - 13:5? 

STAGS? DATA F T?anSvE 3 3E TEE STIFFENS 3 ( U Q0 A L A N a i_ T 5 I 3 ) 

200000 100S3-1 

S 0 0 4 5 3-2 

10253-3 

0, 5 C - 1 

1 0 10000 5 0-2 

I 2 200 . 0 . 2 b 0 0 . 3 0-3 

II 6 11 11 11 3 11 3 11 3 $ F - 1 

1 2 2 a S G-l-1 C0 y3 ATA3ILI TY 1-2 

1 2 3 a s G - 1 - 2 CO w ? A r ABILITY 1-3 

3 2 a a 5 G-l-3 COM? A T AS I L I T Y 3-? 

3 2 S 2 3 G-l-4 C0 m ?ATA3ILITy 3-5 

1 0 S 1-1 

10, Eb ,3 0. . 0?9 3 1-2 MATERIAL PROPERTIES 
1113 K-l-1 FLOOR SHELLS 
1 . 187S 3 K-2-1 

2113 K-l-2 STIFFENS? SHELLS 

1 ,125 3 Y -2-2 

2 0 3 m- i - i 

0, 6. 0. 3. 3 M - 2 A - 1 

1 0 0. -.0033 5 M-5-1 

4 10 3 N - 1 - 1 

a 6 2 a J P-l-1 BOUNDS 1 
0 3 Q- 1 - 1 

0 $ P-l-1 

2 0 3 M-i -2 

0, 6, 3, P. 3 M-2A-? 

1 0 0, - , OP 33 5 M - 5 - 2 
4 10 3 N - 1 - 2 

a 2 2 6 3 P-1-2 BOUNDS 2 
0 S 0-1-2 

0 5 P-1-2 
223 m-i-3 

0, 6, 0. ,7313 3 M - 2 A - 3 

1 3 0, 3. .7313 5 M-3-3 

2 0 PO. 0. 5 m- 5-3 
a 1 0 3 N- 1 -3 

46263 P-1-3 BOUNDS 3 
0 3 Q - 1 - 3 
0 S R - 1 - 3 
2 2 3 M- 1 -4 

0. 6. 3. 3.35P4 3 M -2A -4 

130. 3.35P4 .7313 5 M-3-4 

200. .0625 5 M-5-4 

4 10 3 N - 1 - 4 

43265 P-1-4 BOUNDS 4 
0 3 0-1-4 

0 3 ? - 1 -u 
2 2 3 M- 1 -5 

0. 6. 2.640o 3. 3 m - 2 A - 5 

110. 2. 6406 .7313 I 3-5 
2 0 0. .0625 3 M - 5 - 5 

4 10 3 N - 1 - 5 

4 & 2 3 3 P-1-5 BOUNDS 5 
0 5 0-1-5 
0 3 P-1-5 
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STAG51 DATA Fur? TRANSUE-Sh TEE STIF.-ENER (PLASTIC) 
61 1 100101*3-1 

50041 3-2 
1025 3-3 

0. 5 C - 1 

0 1000 11 * 0-1 

0 2. E-3 20 E~o 5 E- 1 

045 £-2 

1 1 6 1 1 1 1 1 1 3 1 1 3 1 1 3 5 F- 1 
1224$ G- l - 1 COMPATAE ILITY 1-2 
1234$ G-l-2 CGMPATAB ILITY 1-3 
32441 G- 1-3 CGMPATAB ILITY 3-4 

3 2 5 2 $ G- 1-4 C0MPATA3 ILITY 3-5 
13151-1 

10. E6 3 0. .099 $ 1-2 MATERIAL PROPERTIES 

.0041 41 £3 .06 45. E3 1. 47. E3 $ 1-3 PLASTICITY 

1115$ K— 1-1 FLOGR SHELLS 

1 1875 0. 2 $ K-2-1 

2 1 1 5 $ K- 1-2 STIFFENER SHELLS 

1 . 125 0. 2 $ K-2-2 

2 0 $ tt — 1-1 

0. 6. 0. 3. $ M-2A-1 

10 0. 0. 0 1 $ M-5-1 

410 S N— 1—1 

4 6 2 4 * P-1-1 BOUNDS 1 

0 5 Q — 1—1 

0 $ R-l-1 

2 0 $ M- 1-2 

0. 6. 3. 9 $ M-2A-2 

1 0 0 0. 0 1 $ n-S-2 

410 $ N- 1-2 

4 2 2 6 $ P-1-2 BOUNDS 2 
0 $ G-l-2 
0 $ P-1-2 

2 2 $ M- 1-3 

0 6. 0 9375 $ M-2A-3 

1 3 0. 3. 9375 $ 11-3-3 

2 0 90 0. 0 1 $ M-5-3 

410 $ N- 1-3 

46265 P-1-3 BOUNDS 3 
0 5 Q-l-3 
0 5 R-l-3 
2 2 5 M-l-4 

0 6. 3 3. 3594 5 M-2A-4 

1 3 0. 3. 35=54 . 9375 $ rt-3-4 

200 0 015 M-5-4 

410 $ N-l-4 

4 3 2 6 $ P-1-4 30UNDS 4 
0 $ Q-l-4 
•0 $ R - 1-4 

2 2 $ M-l-S 

062. 6406 3 5 M-2A-5 

110 2. 6406 9375 5 M-3-5 

2 0 0 0. 0 1 $ M— 5-5 

410 5 N- 1-5 

4 6 2 3 $ P-1-5 BOUNDS 5 
0 i Q-l-5 
0 t R-l-5 



158 



FLUID TASS DATA PCR TPa;,SvEF8E TEE 5 

0 2^2 150 0 0 
1. 1 52E-3 59 E3 
T T F T 

T F F T 
F F T F 
F F T F 
F T F F 
F T 

TEE. FLU TEE. GEG TEE. COR 
10 

1. 1. 1. 1. 

MESH 

1 4 0 0 0 
17 8 2 

2 4 0 0 0 

2 8 9 3 

3 4 0 0 0 

3 9 10 4 

4 4 0 0 0 
410115 

5 4 0 0 0 

5 11 12 6 

6 4 0 0 0 

7 13 14 8 

7 4 0 0 0 

8 14 15 9 

8 4 0 0 0 

9 15 16 10 
9 4 0 0 0 



10 


16 17 


1 1 


10 


4 0 0 


0 


1 1 


17 18 


12 


1 1 


4 0 0 


0 


13 


19 20 


14 


12 


4 0 0 


0 


14 


20 21 


15 


13 


4 0 0 


0 


15 


21 22 


16 


14 


4 0 0 


0 


16 


22 23 


17 


15 


4 0 0 


0 


17 


23 24 


18 


16 


4 0 0 


0 


L 9 


25 26 


20 


17 


4 0 0 


0 


20 


26 27 


21 


L 8 


4 0 0 


0 


21 


27 23 


22 


19 


4 0 0 


0 


22 


23 29 


23 


20 


4 0 0 


0 


23 


29 30 


24 


21 


4 0 0 


0 


25 


31 32 


26 


22 


4 0 0 


0 


26 


32 23 


27 


22 


4 0 0 


0 


— ~7 
C. / 


32 34 


23 


24 


4 0 0 


0 
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29 


34 35 


29 


25 


4 0 0 


0 


29 


25 36 


30 


29 


4 0 0 


0 


31 


37 39 


32 


27 


4 0 0 


0 


32 


38 39 


33 


29 


4 0 0 


0 


33 


29 40 


24 


29 


4 0 0 


0 


34 


40 41 


25 


30 


4 0 0 


0 


35 


41 42 


36 


31 


4 0 0 


0 


37 


43 44 


38 


32 


4 0 0 


0 


38 


44 45 


29 


33 


4 0 0 


0 


39 


45 46 


40 


34 


4 0 0 


0 


40 


46 47 


41 


35 


4 0 0 


0 


41 


47 48 


42 


38 


4 0 0 


0 


43 


49 50 


44 


37 


4 0 0 


0 


44 


50 51 


45 


38 


4 0 0 


0 


45 


51 52 


46 


39 


4 0 0 


0 


46 


52 53 


47 


40 


4 0 0 


0 


47 


53 54 


48 


41 


4 0 0 


0 


49 


55 56 


50 


42 


4 0 0 


0 


50 


56 57 


51 


43 


4 0 0 


0 


51 


57 53 


C -l 

w — 


44 


4 0 0 


0 


52 


58 59 


53 


45 


4 0 0 


0 


53 


59 60 


54 


46 


4 0 0 


0 


55 


6 1 62 


56 


47 


4 G 0 


0 


56 


62 62 


57 


43 


4 0 0 


0 


57 


o3 64 


58 


49 


4 0 0 


0 


59 


64 o5 


59 


50 


4 0 0 


0 


59 


65 6o 


60 


51 


4 0 0 


0 


6 12 77 67 


52 


4 0 0 


0 


67 


77 7G 


69 


53 


4 0 0 


0 


69 


78 79 


69 


54 


4 0 0 


0 


69 


7 9 80 


70 
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55 

70 

56 

71 

57 

72 

53 

73 

59 

74 

60 

75 

61 

12 

62 

77 

63 

78 

64 

79 

65 

80 

66 

81 

67 

82 

68 

83 

69 

84 

70 

85 

71 

18 

72 

87 

73 

88 

74 

e9 

75 

90 

76 

91 

77 

92 

78 

93 

79 

94 

80 

95 

81 

24 

82 

97 

83 

98 

84 

9 <? 

85 



4 0 0 0 
50 81 71 
4 0 0 0 

o 1 c. / 2 

4 0 0 0 

82 83 73 
4 0 0 0 

83 84 74 
4 0 0 0 

84 85 75 
4 0 0 0 
35 86 76 
4 0 0 0 
18 87 77* 
4 0 0 0 
87 38 78 
4 0 0 0 
38 89 79 
4 0 0 0 

89 90 50 
4 0 0 0 

90 91 81 
4 0 0 0 

91 92 82 
4 0 0 0 

92 93 83 
4 0 0 0 

93 94 84 
4 0 0 0 

94 95 85 
4 0 0 0 

95 96 86 
4 0 0 0 
24 97 87 
4 0 0 0 

97 95 88 
4 0 0 0 

98 99 89 
4 0 0 0 
89 100 90 



4 0 


0 0 




100 
4 0 


101 
0 0 


91 


101 
4 0 


102 
0 0 


92 


102 
4 0 


103 
0 0 


93 


103 
4 0 


104 
0 0 


94 


104 
4 0 


105 
0 0 


95 


105 
4 0 


106 
0 0 


96 


30 
4 0 


107 97 
0 0 


107 
4 0 


1 03 
0 0 


98 


ice 

4 o 


109 
0 0 


99 


l 09 
4 0 


l 10 
0 0 


100 
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101 



100 110 111 
66 4 0 0 O' 

1 0 1 111 112 1 02 
97 4 0 0 0 

102 112 113 103 

aa 4 o o o 

103 113 114 104 
69 4 0 0 0 

104 114 115 105 

90 4 o 0 0 

105 115 llo 106 

91 4000 
30. 36 117 107 



92 4 


0 0 


0 




107 


117 


118 


108 


93 4 


0 0 


0 




108 


1 18 


1 19 


109 


94 4 


0 0 


0 




109 


119 


120 


110 


95 4 


0 0 


0 




1 10 


120 


121 


1 1 1 


96 4 


0 0 


0 




1 1 1 


121 


122 


1 12 


97 4 


0 0 


0 




1 12 


122 


123 


113 


98 4 


0 0 


0 




1 13 


123 


124 


114 


99 4 


0 0 


0 




114 


124 


125 


115 


100 


4 0 


0 0 




115 


125 


126 


116 


101 


4 0 


0 0 




36 42 127 1 : 


L 7 


LG2 


4 0 


0 0 




1 17 


127 


128 


1 18 


103 


4 0 


0 0 




1 18 


129 


129 


1 19 


104 


4 0 


0 0 




119 


129 


130 


120 


105 


4 0 


0 0 




120 


130 


131 


121 


106 


4 0 


0 0 




121 


131 


132 


122 


107 


4 0 


0 0 




122 


132 


133 


123 


108 


4 0 


0 0 




123 


133 


134 


124 


109 


4 0 


0 0 




124 


134 


135 


125 


110 


4 0 


0 0 




125 


135 


1 3o 


126 


1 1 1 


4 0 


0 0 




42 ae 137 i: 


27 


1 12 


4 0 


0 0 




127 


137 


133 


129 


113 


4 0 


0 0 




128 


1 3 8 


139 


129 


1 14 


4 0 


0 0 




129 


139 


140 


130 


1 15 


4 0 


0 0 




130 


140 


14 t 


131 



162 






1 15 


4 


O' 


0 


0 




121 


1 


4 1 


1 


42 


1 22 


117 


*>♦ 


0 


0 


0 




132 


1 


- n 

*“ L CL 


1 


43 


123 


i ia 


4 


0 


0 


0 




132 


143 




44 


124 


119 


4 


0 


0 


0 




134 


144 


145 


125 


120 


4 


0 


0 


0 




135 


145 


146 


126 


121 


4 


0 


0 


0 




48 : 


54 


147 


127 


122 


4 


0 


0 


0 




137 


147 


148 


128 


123 


4 


0 


0 


0 




138 


148 


149 


129 


124 


4 


0 


0 


0 




139 


149 


150 


140 


125 


4 


0 


0 


0 




140 


150 


151 


141 


125 


4 


0 


0 


0 




141 


i 


51 


1 


c n 

w C_ 


142 


127 


4 


0 


0 


0 




142 


152 


153 


142 


128 


4 


0 


0 


0 




143 


153 


154 


144 


129 


4 


0 


0 


0 




144 


154 


155 


145 


130 


4 


0 


0 


0 




145 


155 


156 


146 


131 


4 


0 


0 


0 




54 60 


157 


147 


132 


4 


0 


0 


0 




147 


157 


1 


58 


140 


133 


4 


0 


0 


0 




140 


1 


53 


1 59 


1 49 


134 


4 


0 


0 


0 




1 49 


159 


160 


150 


125 


4 


0 


0 


0 




150 


160 


161 


151 


135 


4 


0 


0 


0 




151 


lol 


162 


152 


137 


4 


0 


0 


0 




152 


162 


163 


152 


128 


4 


0 


0 


0 




153 


163 


164 


154 


139 


4 


0 


0 


0 




154 


1- 


64 


165 


1 55 


140 


4 


0 


0 


0 




155 


165 


l 6c 


1 56 


141 


4 


0 


0 


0 




60 86 


1 A 


1 7 


1 57 


142 


4 


0 


0 


0 




157 


167 


168 


153 


143 


4 


0 


0 


0 




15c 


1 63 


169 


1 59 


1 44 


4 


0 


0 


0 




1 59 


l' 




i ■ 


70 


1 60 


145 


4 


0 


0 


0 




160 


1 70 


17 1 


let 


1 4 £ 


4 


0 


0 


0 





163 



1 6 1 


1 7 1 


172 


162 


147 


4 0 


0 0 




lo2 


172 


1 72 


1 o3 


142 


4 0 


0 0 




162 


172 


174 


164 


149 


4 0 


0 0 




164 


174 


175 


165 


150 


4 0 


0 0 




165 


175 


176 


166 



MASS 

EIGEN 

EXIT 



164 



lL iL ll. 



AlJOMEN TF - J ' -* ~ P I f D 4 ' 4 ^ uP i « Ar-iS'-' ‘r: E T n n bT i - - i'cr 

i ci M a- E i ci. " wJ i z. ‘ j. ; ju - i i- r. ~ w r. 
r 
r 
F 

0 

242 1452 6 3 
1 

0 1 1 50 1 
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TEETi.T lNP 



L9-QC7-19E4 15 

TIME IN7EGPA7CP DATA PGP TP ANE^EPEE TEE STIFFENS 
TEE. P-E TEE. PCS 
TEE P 5 7 

0 G i 05. 7E 

0. 2 S4 9.375 

T F F F 

F 

101 

1 

2E-4 

3125 0 1. 56E-3 
1 

0 20 E-6 
2. E-3 

1 20 

0 0 0 0 
F F 
F 
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O WJ 111 



r-crG^ i l U? , I 2G--jC r - 1 A 12 OA 

\jS a P 0 ji T — p R CC G 2 5 OF Data PGP 7RAr,5VGFSG TEE STIFFENED 
1 

l 

E PRE 
E FOE 
F T F A F 
F F F F 
F F T 

0 G 2 

3 1 o l 

1 67 7o 1 
0 0 2 

2 16 1 

3 ti 7 76 1 
0 2 
1 c, 1 
ol *9 1 
F 

1 0 

F F F F 
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TESS T- 


02 1 


1 I'Ju « 


L PL 


UV i 


*111 


0 








0 0 0 


0 








7 22 


^ r > 

— 


2 3 






0 

- 7 22 


-2 


. 22 






l 0 0 


0 1 


100 


0 0 




7 i 2 


22. 


2 3 


. 3 




2 1 0 


0 6 


1 00 


0 1 


0 0 


. 7 22 . 


22. 


23 . 


10. 


10 


2 1 0 


0 6 


100 


0 2 


0 0 


7 22 


23 


23 


10. 


10 


2 1 0 


0 * 


too 


0 2 


0 0 


. 7 22 . 


22. 


23 


10. 


10 


2 2 0 


0 6 


100 


0 1 


0 0 


7 23 


23 


22 


10 


10 


2 2 0 


0 6 


100 


0 2 


0 0 


7 22 


73 


2 2 


10. 


10 


2 2 0 


0 0 


100 


0 3 


0 0 


. 7 22 . 


2 ■“> 


23 . 


10. 


10 
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t 



RECSTAGS2. I'JPMO 



9 - 5 EP- 195 J 



10:53 



ST A OS 2 Di'i r DP LONGITUDINAL RE Z T A NGUL A R ST: rr £NER 

500000 t 0035— i 

3 0 0 2 3 3-2 

1 0 1 S 8-3 

0. S C-l 

l 0 10000 5 0-2 

12 200. 0. 2300. S 0-3 

6863355F-1 

l 3 2 l 5 G - 1 - L COH° A T A3 I L I T Y 1-2 

l 3 3 1 5 G-l-2 COV= AT ABILITY 1-3 

1031-1 

lO.OEb .3 0. .005 5 1-2 MATERIAL PROPERTIES 

1 l 1 $ K- 1 

t .25 3 K-2 

2 0 5 M- l - l 

0. b. 0. 0. 5 - 2 A - l 

1 0 0 . -0.125 3 *-5-1 
a 1 0 5 N- 1 -l 

2 2 6 a 3 P - l-l BOUNDS 1 
0 3 Q - l - I 

0 S R-l -1 

2 0 S H-i -2 

b. 12. 0. <?. 5 M-2A-2 

1 0 0 . - 0.125 S M- 5-2 
a i 0 S N - 1 - 2 

b 2 2 U S P-1-2 BOUNDS 2 
0 5 0 - 1-2 
0 S R-l -2 

2 2 3 *- 1 -3 

0 . 1.25 0 . 9. 3 M - 2 a - 3 

3 l 6. 0. 1.25 3 M-3-J 

1 0 90 . 0 . S vi- 5-3 
a t 0 3 N - 1 - 3 

b 2 3 a 3 P-1-3 BOUNDS 3 
0 3 0-1-3 
0 $ R-l -3 



( '•'COAL 



r 2 I 5 j 
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ETaGSi DATA FGP UGNO l r jZ :NA W R EG TanGULaF 
c 1 1 l GO 1 G- 1 5 3-1 

3 G G 2 5 5-2 
1015 B-3 
G. 5 C - l 

0 1C0G 1 1 5 D-l 

0 2. £-3 20 E-6 s E- 1 
Gas £-2 
6 0 6 3 3 3 5 F— 1 
13215 G-l-l CGMPATABILITY 1- 
13 3 15 G-l-2 COMPARABILITY 1- 

1 3 1 s I- 1 - - -- - -- 

10. 0E6 .3 0. .098 5 1-2 MATERIAL PROFESTI 

. 004 1 41. E3 .06 45. E3 1. 47. E3 5 1-3 PLAS 

1115 5 K-l .*-•• 

1 . 25 0. 2 5 K-2 “ * ' ' 

2 0 5 M-l-1 • • : 

0. 6. 0. 9. 5 M-2A-1 

10 0. 0. 0 15 M-5-1 

410 5 N— 1—1 

22645 P-1-1 SOUNDS 1 
0 5 G-l-l . 

0 5 R-l-1 
205 M- 1-2 
6 12 0. 9. 5 M-2A-2 

10 0. 0. 0 15 M- 5-2 

410 5 N-l-2 

62245 P-1-2 3GUNDS 2 
0 £ 0 - 1 - 2 

0 5 R- 1-2 
225 M-l-3 

0. 1. 375 0 9 5 M-2A-3 

~3 1 6. 0. 1 3~5 * M — 3 — 3 

1 3 90 0. 0 1 5 M-5-3 

1 0.0 9 00 1 

•: 0 0 45. 0005 

'•j 'J , si ■ <=. J 0 Vj 1 D 

410 5 .N-1-2 

- 3 2 -*■ 5 P-1-3 3GUND5 3 
■ 5 0-1-3 
0 5 P.-1-3 
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U M 



R EC T ANGULAR 



FLU I D M A 3 S D A T A F u 9 LONGITUDINAL 

0 104 70 0 0 

1. 152E-3 59 0E3 
T T F 7 
T F T F 
F F T F 
F F T F 
F T F F 
F T 

REG. FLU REC.QEO REC.CER 
10 

i. n. 

o. o. -1. o. 

MESH - .. >. v 

1 4 0 0 0 

1 9 10 2 

2 4 0 0 0 
210113 

3 4 0 0 0 

3 11 12 4 

4 4 0 0 0 

4 12 13 5 ill; 

5 4 0 0 0 

5 13 14 6 

6 4 0 0 0 

6 14 15 7 

7 4 0 0 * 

7 15 163 

8 4 0 0 0 

9 17 13 10 

9 4 o 0 C 

10 18 1* 

10 4 0 0 

11 19 20 
114 0 0 

12 20 21 

12 4 o 0 

13 21 22 

13 4 0 0 

14 22 23 

14 4 0 0 

15 23 24 

15 4 0 0 
17 25 26 

16 4 0 0 
13 26 27 

17 4 0 0 

19 27 23 
13 4Q0 

20 23 29 

19 4 Q o 

21 29 30 

20 4 0 0 

22 20 31 

21 4 0 0 

23 31 32 

22 a o 0 

25 33 34 

23 4 0 0 

26 34 35 



1 1 
0 
12 
0 

13 
0 

14 
0 

15 
0 

16 
0 
18 
0 

19 
0 

20 
0 
21 
0 

22 

0 

23 

0 

-1 

C . ~ r 

0 

26 

0 

27 



S7IFFNE9 
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24 


o 

o 


0 


27 


3 5 36 


23 


25 


4 o 0 


0 


29 


36 37 


29 


26 


4 0 0 


0 


29 


37 38 


30 


27 


4 0 0 


0 


30 


38 39 


31 


28 


4 0 0 


0 


31 


39 40 


32 


29 


4 0 0 


0 


33 


41 42 


34 


30 


4 0 0 


0 


34 


42 43 


35 


• 31 


4 0 0 


0 


35 


43 44 


36 


32 


4 0 0 


0 


36 


44 45 


37 


33 


4 0 0 


0 


37 


45 46 


39 


- 34 


4 0 0 


0 


28 


46 47 


39 


35 


4 0 0 


0 


39 


47 48 


40 


36 


4 0 0 


0 


41 


49 50 


42 


37 


4 0 0 


Q> 


42 


50 51 


42 


38 


4 0 0 


0 


43 


5 1 52 


44 


39 


4Q0 


0 


44 


52 53 


45 


40 


4 0 0 


0 


45 


53 54 


46 


41 


4 0 0 


0 


46 


54 55 


47 


42 


4 0 0 


0 


47 


55 56 


43 


43 


4 0 0 


0 


49 


57 58 


50 


44 


4 0 0 


0 


50 


58 59 


51 


45 


4 0 0 


0 


51 


59 60 


52 


46 


4 0 0 


0 


52 


60 61 


53 


47 


4 0 0 


0 


53 


ot 62 


54 


48 


4 0 0 


0 


54 


62 63 


c = 


49 


4 0 0 


0 


55 


o3 64 


56 


50 


4 0 0 


0 


57 


65 66 


58 


51 


4 0 0 


0 


58 


66 6 7 


59 


52 


4 0 0 


0 


5 9 


67 o6 


oO 


3. j 


‘*0 0 


0 


60 


68 69 


6 l 


5 ; 


4 0 0 


0 
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61 69 70 
5 5 4Q0 

62 70 71 

56 4 0 0 
o2 71 72 

57 4 0 0 
65 73 74 
56 4 o 0 
6o 74 75 

59 4 0 0 

67 75 76 

60 4 0 0 

68 "76 77 

61 4 0 0 

69 77 78 

62 4 0 o 

70 79 79 

63 4 o 0 

71 79 80 

64 4 o o 

73 81 82 

65 4 o 0 

74 82 83 

66 4 0 0 

75 e3 84 

67 4 0 0 

76 84 85 

68 4 0 0 

77 85 86 

69 4 0 0 

78 86 87 

70 4 0 0 

79 87 88 
MASS 
EIGEN 
EXIT 



62 

0 

63 

0 

6 4 

0 

co 

0 

67 

0 

68 

0 

69 

0 

70 

0 

71 

0 

72 

0 

74 

0 

75 

0 

76 

0 

U 

0 

73 

0 

79 

0 

80 



IJL tL I- 



AUG 

°EC 



0 

IGA 

1 

0 1 



’EN TED flA T^ 3 I 
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